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FOREWORD 

The  attractive  prospect  of  estimating  wind  profiles  directly  from 
satellite  radiance  data  is  motivated  by  three  factors: 

•  the  abundance  of  satellite  data  over  areas  such  as  the 
oceans  which  lack  conventional  data  sources 

•  the  cost-effective  advantages  of  obtaining  representative 
winds  in  data-sparse  regions  without  additional  instru¬ 
mentation,  and 

•  the  avoidance  of  errors  associated  with  an  indirect 
approach  based  on  collections  of  temperature  profiles 
obtained  from  satellite  radiance  data  by  solving  the  inver¬ 
sion  problem. 

Several  attempts  to  develop  a  direct  approach  have  been  discussed 
in  the  open  literature  but  have  yielded  disappointing  results.  The  most 
important  of  these  efforts  are  reviewed  in  the  proposal  which  led  to  this 
project  (Ref.  2)  which  is  funded  by  The  Satellite  Meteorology  Branch  of 
the  Air  Force  Geophysics  Laboratory.  As  discussed  in  Appendi.x  G,  an 
interim  technical  report,  this  project  extends  previous  work  by  minimizing 
a  specified  error  criterion  to  obtain  statistic.  '  >  optimal  thermal  wind  esti¬ 
mates.  The  framework  established  for  incorp  -  uing  statistical  error  mod¬ 
els  of  relevant  data  sources  and  algorithm  components  is  also  a  significant 
contribution  of  the  project. 
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EXECLTIVE  SUMMARY 


This  repon  summarizes  a  research  project  conducted  by  TASC  during  GF>'88.  with  ap¬ 
proximately  one  person-year  level  of  effort.  I'he  central  idea  behind  the  project  was  to  use  sta¬ 
tistical  models  of  error  sources  and  an  optimal  estimation  theory  approach  to  estimate  thermal 
aim!  profiles  directiv  from  multi-spectral  satellite  sourider  radiance  data.  (VVe  refer  to  the  proce¬ 
dure  as  the  optimal  estimation  algorithm  in  this  repon.)  Thermal  wind  profiles  contain  esti¬ 
mates  of  the  vector  difference  between  the  geostrophic  wind  at  each  of  several  altitudes  and  that 
at  a  single  reference  altitude.  The  reference  altitude  would  typically  be  at  the  lower  end  of  the 
podile,  e  g.,  that  of  a  surface  wind  measurement  by  a  ship-based  instrument.  The  project  was 
motivated  by  a  need  for  im[iroveil  wind  measurements  in  data  sparse  regions  such  as  cloud-free 
aretis  over  the  ocean  If  successful Iv  Implemented,  the  resulting  algorithm  would  provide  a  com¬ 
plement  to  cloud-drift  wind  estimation  procedures.  Figure  ES-1  contains  a  schematic  represen¬ 
tation  of  how  an  operational  (Optimal  estimation  algorithm  might  function. 

.Section  2  of  this  Final  Report  tlemonstrates  clearly  that  the  algorithm  generates  reason¬ 
able  thermal  wind  profile  estimates  using  only  approximate  first-guess  information  typical  of 
what  would  exist  in  an  open  ocean  region  Figure  FS-2  presents  results  for  two  initial  data  con¬ 
ditions  along  with  raw  insonde-measured  values  of  total  wind  shear  near  the  estimate  location. 
The  two  initial  data  conditions  are:  1)  using  a  temperature  derivative  profile  derived  from  an 
analysis  field  based  on  rawinsonde  data  and  2)  using  a  crude  approximation  to  the  analyzed 
profile  typical  of  what  a  forecast  or  analysis  field  might  provide  at  an  open-ocean  location.  For 
each  case,  initial  profiles  (dashed  curves)  and  final  estimates  (solid  curves)  are  shown.  The  im¬ 
portant  point  demonstrated  by  this  figure  is  that  even  though  the  initial  data  provided  much  dif¬ 
ferent  first-guess  wind  profiles  (calculated  from  the  temperature  derivative  profiles  via  the 
thermal  wind  equation),  the  algorithm  produced  very  similar  final  estimates  from  tlte  satellite 
radiance  data. 

The  total  shear  data  plotted  in  Fig.  ES-2  are  rawinsonde-measured  east  wind  differences 
between  each  pressure  level  and  850  mb.  .Although  the  shape  of  the  total  shear  profile  is  similar 
to  that  of  the  other  curves,  it  is  clear  that  a  significant  non-geostrophic  wind  component  existed 
in  the  region.  1  his  is  clearly  a  limitation  of  any  wind  estimation  approach  relying  solely  on 
radiance  measurements. 
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Figure  ES-l  Oveniew  of  Optimal  Thermal  Wind  Profile  Estimation 

Analyses  described  in  this  report  ^ho\v  that  significant  data  sensitivities  may  exist  which 
could  make  superiority  of  the  optimal  estimation  approach  over  alternate  approaches  unclear. 
For  example,  a  primary  alternate  approach  is  to  develop  (i.e..  estimate)  a  3-D  temperature  field 
from  a  grid  of  satellite  soundings  by  "temperature  inversion  "  applied  to  each  sounding.  Then 
the  themial  wind  equation  can  be  used  to  estimate  the  thermal  wind  profile  from  the  layer-aver¬ 
aged.  horizontal  gradients  of  the  temperature  field  The  computational  burden  of  this  indirect 
approach  is  significantly  larger  than  that  of  direct  optimal  estimation,  due  primarily  to  the  re¬ 
quirement  for  N’  temperature  inversions  (for  an  N  \  N  grid  of  radiance  data),  and  is  clearly  sus¬ 
ceptible  to  inversion  errors  through  the  numerical  derivatives  required  to  define  horizontal 
temperature  gradients. 

However,  the  optimal  estimation  approach  also  requires  gradient  data,  in  this  case  the 
horizontal  gradients  of  the  radiances  at  each  sounder  wavelength.  Experience  with  the  sounder 
data  used  in  the  demonstration  showed  that  these  gradient  estimates  can  be  more  sensitive  to 
the  size  and  location  of  the  data  grid  than  is  predicted  by  error  models  derived  from  straight- 


ES-2 


G-15592 

5-l9-«9 


Figure  KS-2  T  hermal  Wind  Profile  Estimates  by  the  Optimal 
I'.stimation  Algorithm 


forward  radiance  data  models.  Recommendations  in  Section  3.2  focus  on  examining  further, 
and  reducing  sensitivity  to.  the  data  variability  through  additional  error  model  development.  If 
these  reductions  in  sensitivit)  are  achieved,  then  the  optimal  estimation  algorithm  will  provide  a 
superior  approach  to  thermal  wind  profile  estimation  relative  to  both  the  indirect  approach  de- 
serihed  abo\c  and  to  alternate  direct  approaches  discussed  in  Refs.  1  and  2. 
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1. 


INTRODUCTION 


This  Final  Repon  of  rAS(”s  Wind  F’rofile  Estimation  project  provides  a  summary  of 
theoretical  and  nunierical  results  obtained  ^ince  delivery  of  the  Interim  Repon  (the  complete  in¬ 
terim  report  is  included  in  this  volume  as  Appendix  G).  At  an  interim  project  review  meeting 
held  at  AFCiL  in  April  19.S,S,  much  discus^ion  focused  on  the  imponance  of  determining  the  sen¬ 
sitivity  of  estimates  to  various  model  elements.  .Accordingly,  project  emphasis  was  shifted  from 
model  development  to  sensitivitv  analysis  to  identify  those  model  components  most  strongly 
influencing  estimation  accuracy.  In  this  way.  future  effons  can  focus  on  the  dominant  eiTor 
contributors. 

Most  activity  since  the  Interim  Report  and  meeting  has  been  in  three  areas:  evaluating 
the  new  wind  profile  estimation  algorithm  iiNing  the  sample  data  set;  conducting  sensitivity 
analyses  of  the  algorithms;  and  comparing  the  new  algorithm  with  Ohring's  method  (Ref.  1). 
both  as  to  implementation  and  numerical  ptxtperiies  and  to  parameter  sensitivity.  Results  in  all 
of  these  areas  are  cttntained  in  this  technical  report. 

Section  2  of  the  repon  contains  four  subsections  which  present  the  most  imponant  re¬ 
sults  of  the  project.  T  he  first  subsection  describes  calculation  of  radiance  gradient  vectors  from 
the  sample  data  set.  Issues  of  gradient  variability  with  calculation  grid  size  and  location  are  de¬ 
scribed  in  detail  (see  also  Appendices  .A.  B,  and  D).  Results  related  directly  to  the  new  optimal 
estimation  algorithm  are  presented  in  the  second  subsection  of  Chapter  2  and  Appendices  C 
and  E.  These  include  an  interesting  and  imponant  theoretical  result  which  shows  that  the  algo¬ 
rithm.  as  implemented,  scales  the  first-guess  wind  profile  so  as  to  best  fit  the  radiance  field  data. 
Numerical  results  (i.e.,  thermal  wind  profile  estimates)  resulting  from  several  initial  tempera¬ 
ture  gradient  profile  assuniptions  are  presented  in  this  subsection.  The  third  subsection  presents 
results  of  apply  ing  O’  •■.j,'s  -.nethod  (Ref.  1)  to  the  same  data  set.  Finally,  the  numerous  sensi¬ 
tivity  results  are  Ccht  .d  in  the  fourth  subsection  of  Chapter  2  and  Appendix  F. 

Delivery  of  thi'-  un  constitutes  completion  of  research  effons  for  this  program.  .A 
summaiT  of  the  co 'I'.'lusions  which  can  be  drawn  is  contained  in  Chapter  3.  While  all  of  the 
tasks  outlined  in  the  contract  Statement  of  Work  have  been  completed,  additional  effon  will  be 
required  to  fully  implement  a  temperature  gradient  profile  error  model  for  the  algorithm.  That 


ett'on  and  other  task>  leading  toward  an  operational  algoritlim  .11  c  doeiissed  briefly  in  the  fintil 
subsection.  Those  tasks  should  letid  to  a  practical,  operational  algoi  ithni  for  producing  thermal 
wind  profiles  directly  from  satellite  radiance  data. 
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RESULTS 


2.1  R\DIA.\CE  GR.\I)IE.M  MXTORS 

High  Resolution  Intiared  Radiation  Sounder  2  (HIRS-2)  data  selected  for  analysis  were 
from  the  11x11  grid  shoun  in  the  box  in  f  ig  2.1-1.  Contour  lines  in  the  figure  are  radiance  in¬ 
tensities  for  the  visible  channel  (channel  20).  the  lx)x  is  located  in  a  region  of  relatively  low  in¬ 
tensity  from  this  chtinnel,  ini.licating  ;i  likel\  clouii-fiee  tiiea.  (This  is  consistent  with  AVHRR 
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Figure  2.1-1  Visible  Channel  Intensity  Contours  and  Data  Region 


images  ot  the  recion.)  Fieure  2  1-2  shows  the  grid  location  uiihiii  a  geopolitical  map  ot  the 
western  Mediterranean  region 

In  Appendix  A.  plots  ot  the  data  from  the  six  channels  nsc-d  in  this  demonstration  are 
presented.  Table  2.1-1  identifies  the  channels  and  their  center  wavelengths.  There  are  two  plots 
in  Appendix  .A  for  each  channel,  each  containing  raditince  values  aU>ng  six  scan  lines  (the  cen¬ 
ter  of  the  1  1  scan  lines  is  in  both  plots).  Scan  lines  and  track  axis  distances  are  identified  by  lo¬ 
cal  track-scan  coordinates  (in  km)  in  which  (0.  0)  is  at  the  northwest  corner  of  the  region  and 
(4(10.  400)  at  the  southeast  corner.  (Note  that  the  irack-scan  axes  are  rotated  at  -  20°  rela¬ 
tive  to  nonh-east  coordinate'.! 

Estimates  of  the  radiance  gradients  for  all  ch.tnnels  were  calciikited  from  local  planar 
fits  to  the  data  surface  within  the  region  of  interest.  .Appendix  I)  presents  details  and  an  error 
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Figure  2.1-2  Geographical  Location  of  Data  Region 


Table  2.1-1  Suniinarv  of  HIR.S-2  Data  used  iji  the  Demonstration 
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inojfl  for  the  eraJieiit  e^Ii^laIi^)n  procciluix  ( iiadient  e-'iimaies  depend  criiicnlly  on  the  region 
over  which  the  planar  fit  is  made  Region  selection  was  based  on  observed  variability  ot  the  esti¬ 
mated  gradients  with  region  size  and  location. 

Figure  2. 1 -a  contains  the  rttdiance  giadient  vectors  computed  at  tlie  grid  center-point  ns- 
ing  5  grid  sizes  for  channels  3  ,ind  1.''.  Cleaiiy  the  apparent  slope  at  the  center  depends  verv 
much  on  tlie  scale  over  which  the  slope  is  evaluated.  Fixamination  of  vector  plots  similar  to 
Fig.  2  1-3  for  all  si\  channels  reveals  the  following: 

•  I  sing  the  3  \  3  grid  vields  gradients  much  different  than  those  from  the 
.s  \  and  larger  grids 

•  L  sing  3  \  and  "  \  g'  ivls  yields  similar  gradients 

•  Using  9  \  d  and  11  x  11  grid  vields  similar  gradients 

•  The  y  .\  y  gradients  are  usually  quite  different  than  the  7  \  7  gradients. 

These  observations  indicate  that  the  d  x  d  (and  1 1  .\  II)  grids  centered  at  the  center-point 
of  the  data  region  contain  data  influenced  by  clouds.  Presence  of  those  clouds  vvas  also 
indicated  by  the  relatively  higher  channel  20  returns  near  the  southeast  and  nonheast  comers  of 
the  data  region  (see  f'ig.  2.1-1). 

To  quantify  the  variability  of  the  radumce  gradients  with  position  of  a  7  x  7  calculation 
grid,  mean  values  and  sttindard  deviations  tor  each  channel  as  the  7x7  grid  is  moved  over  the 
25  possible  locations  within  the  11x11  data  region  are  shown  in  Table  2.1-2.  along  with  the 
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corresponding  model-based  standard  error  for  the  gradient  component  estimates.  Sample  stan¬ 
dard  deviations  in  Table  2  1-2  iiKlicate  that  there  is  n/tt  variability  for  channels  4.  5.  hand 
15  as  the  center-point  of  the  7  x  arrav  is  moved  about  the  data  region.  Some  of  that  variability 
is  certainly  due  to  the  cloud  elfects  alreadv  noted  near  the  nonheast  and  southeast  corners  of 
the  region.  The  sample  standai  d  deviations  are  substantially  larger  than  the  model-based  stan¬ 
dard  deviations  stiggesting  that  [■'rcsence  of  the  clouds  causes  non-negligible  "model-error  "  in 
the  planar  fit  approach  to  radiance  gradient  estimations  (see  the  discussion  in  .-Xppendi.x  D.2). 

After  some  experimentation,  and  consistent  with  the  idea  of  using  as  large  an  area  as 
possible  on  which  to  base  the  gradients  (whicli  are  usei-l  by  the  s\ noptic-scale  thennal  wind 
equation),  it  was  decided  it)  use  graelieiiis  computed  using  '  \  grids  L’nless  stated  otherwise, 
all  estimation  results  presented  in  the  rest  of  this  section  are  based  on  the  radiance  gradients  es¬ 
timated  at  the  center  t)f  the  11x11  data  region  using  a  "  x  "  region  for  the  planar  fit.  These  gra¬ 
dients  tire  tiibulated  for  all  six  channels  m  1  alile  2  1-3  along  w ith  the  standard  dev iation  of  the 
estimation  error  for  each  channel  baseel  on  the  error  model  presented  in  .Appendix  D  (see 
Eq.  D.2-8J. 

Appendix  B  contains  plots  showing  the  normalized  grtidiont  vectors  based  on  7  x  pla¬ 
nar  fits  centered  ;u  each  of  d's  grid  points  for  each  t)f  the  six  channels.  .Although  there  is  still 
some  variability  in  the  gixidient  direction  <inu)ng  the  computed  values,  they  do  not  vaiy  wildly 
from  one  grid  point  to  tinother. 


Table  2.1-3  Kstiinated  Radiance  (iiadients  at  Center  of  Data 
Collection  Region  (7x7  Planar  Fit)''' 


(1)  t  nits  -lie  (niVV  III)  .  (m-  i.m'’  ^ri 


OPTIMAL  ESTIMATION  METHOD 


■>  T 

2.2.1  Optimal  Estimation  .Algorithm  Interpretation 

•Althoiigli  no  t'undanientally  new  equations  liave  been  developed  which  require  inodilica- 
tion  to  the  description  of  the  optimal  estimation  algorithm  found  in  Section  3  of  the  Interim  Re¬ 
port  (.Appendi.\  G).  implementation  ot  the  algorithm  has  led  to  additional  insight  regarding  its 
interpretation. 

Consider  estimation  of  the  east  component  of  the  thermal  u  ind  shear  (relative  to  a  fi.xed 
lowest  le\el.  Zi)  at  each  ot  n  levels; 

(2.2-1) 

(The  north  compv-'nent  can  he  treated  in  a  similar  lashion.)  .As  shoun  in  .Appendi.x  (i.  this  sector 
of  estimates  can  be  computed  using 

u  ---  -(AVI)  (  g  (2.2-2) 

where  g  is  a  sector  of  partial  derisatises  of  the  radiances  w  ith  respect  to  the  y  coordinate  and  R 
and  f  are  the  gas  and  Coriolis  constants  respectisely.  Hie  optimtil  set  of  sseighting  constants  is 

C  =  (I  (f'  P  '  (2.2-3) 

where 

P  =  I'  +  0  o'  (2  2-4) 

0  =  f''  K(c)  ''''  (Ode  (2-2-5) 

J  n  <>y 

and 

rt,  -  [  —  (C)  dC  i  =  I . n  (--2-6) 

J.'  '*y 

In  Eq.  2.2-5.  z  =  In  (pro/P)  where  p  is  the  pressure  level  and  Pri,f  is  an  arbitrary  reference  pres¬ 
sure.  Zf  is  the  z  value  at  the  (appro.ximate)  top  of  the  atmosphere  (zj  corresponded  to  the 
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10  mb  level  in  all  compiitaiionsj  and  K(0  i''  the  vector  of  kernel  functions  defined  in  Appen¬ 
dix  Ci 


By  rearranging  the  tactors  in  the  optimal  estimator,  Eq.  2.2-2,  the  vector  of  estimates  can 
be  uritten; 

u  -  Ig'  B  '  0)  (-  (R/fj  <±)  (2.2-7) 

But  the  second  factor  is  the  estimaie  one  'Aould  compute  dii'cctly  from  the  assumed  temperature 
graiiicnt  profile  using  the  thermal  uiikI  equatmn. 

u’  =  -(A’/f)  (i  (2.2-S) 

So,  the  optimal  uind  pridile  estimate  Ciin  be  written 

u  d  u"  (2.2-9) 

where  d  is  the  scalar  constant  found  in  the  first  factor  of  (:q.  2.2-''  which  does  not  depend  on  the 
leeels  ttt  which  compiments  of  u  tire  defined. 

Now.  examination  of  Etp  2.2-0  shows  cletirh  thtit  the  optimal  estimate,  bused  on  all  of  the 
satellite  radiance  data,  is  simpl\  a  scaled  version  cif  the  initial  thermal  wind  estimate.  .-Xs  will  be 
brought  out  in  the  discussion  of  results,  tlie  tissumed  temperature  gradient  profile  establishes 
the  shape  of  the  vertictil  profile;  then  the  optimtil  estimation  procedure  uses  the  radiance  gradi¬ 
ent  data  to  scale  that  profile  (b\  d)  so  as  to  best  explain  those  observed  gradients. 

Equation  2.2-2  is  based  on  the  error  model  for  the  temperature  gradient  estimates  but 
not  on  the  error  in  the  tissumed  temperature  gradient  profile.  Therefore  the  above  interpretation 
does  not  esttiblish  an  absv)lute  limit  on  the  usefulness  of  the  satellite  data  for  this  application. 
Once  a  temperature  gradient  [irofile  error  model  has  been  incorporated,  provided  that  those  er¬ 
rors  (i.e,.  the  variances)  change  with  altitude,  some  reshaping  of  the  thermal  wind-derived  pro¬ 
file  can  be  expected. 

2.2.2  Temperature  Data  and  .Vpproximalions 

I  he  optimal  estimation  algorithm.  Eq.  2.2-2.  is  based  on  assumed,  measured  or  esti¬ 
mated  temperature  and  temperature  gradient  profiles.  The  Gibraltar  sounding  taken  one  hour 
before  the  satellite  radiance  measurements  were  acquired  vvas  used  for  this  demonstration. 
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I ProceJuie^  uliich  could  be  u>ed  in  an  operational  implementation  are  discussed  briefly  in  Sec¬ 
tion,''.!  The  Intei'im  Repoit,  .Appendix  G.  described  the  satellite  radiance  data  set,  the 
rauindsone  data,  and  the  s\noptic  situation  on  which  this  demonstration  is  based  in  some  detail 
Figure  2  2-1  contains  the  temperature  profile  generated  from  the  rawinsonde  sounding;  the 
Planck  functions  used  in  calculating  the  kernel  functions  were  computed  using  this  profile.  Plots 
of  the  kernel  functions  for  each  channel  are  contained  in  .Appendix  C.  Those  profiles,  and  sev¬ 
eral  others  shown  in  this  report,  are  plotted  against  the  z  Vctriable  defined  in  the  previous  subsec¬ 
tion  with  p.,;  chosen  (arbitrarilv  l  to  be  101.' mb. 

The  areal  co\erage  of  a  temperature  analysis  on  mx  isobaric  surfaces  which  was  used  to 
compute  temperature  gi'adients  is  illustrated  in  t  ie  2  2-2.  The  figure  contains  the  S5()  mbNMC 


r \2'H  ? 
12  ir  " 


PRESSURE 

(MB) 


TEMPERATURE  (DEG  C) 


Figure  2.2-1  Gibraltar  Temperature  Profile 
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Fi^iure  2.2-2  N\IC,  Tciiipcrature  Anaivsis  (S5()  mb)  (deg  C) 

temperature  field  for  1200  L  I’C  4  March  10, S2  interpolated  to  a  1  deg  by  1  deg  grid  ( —  100  kni) 
and  contoured  in  2  deg  (  increments  The  horizontal  temperature  gradient  on  each  surface  was 
computed  at  every  grid  point,  then  values  within  a  5.x5  subgrid  centered  on  Gibralter  were  aver¬ 
aged.  N  ielding  representative  gradients  at  .S.sO,  700.  500,  500,  200  and  100  mb.  (.A  5  x  3  subgrid 
could  have  been  selected,  since  its  area  would  approximate  the  size  of  the  area  over  which  radi¬ 
ance  gradients  were  computed.  However,  it  was  decided  that  the  5  x  5  subgrid  would  provide 
more  information  about  the  atmosphere  around  and  above  Gibraltar.)  The  results  of  the  gradi¬ 
ent  computation  are  shown  in  Table  2.2-1. 


Table  2.2-1  Temperature  (Jradients  Based  on  a  5  x  5  Sul)^rid  (deg  C/ni) 
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The;  temperature  ^radlcnt^  at  m.x  levels  were  iiUerpolateJ  tiiul  extrapokited  to  produce 
the  profiles  slioun  in  Fig.  2.2-.'.  Note  well  tlie  sienifiviinily  larger  magnitude  of  the  north 
component  which  will  be  reflected  throughout  the  report  h\  the  emphasis  on  the  east  component 
of  the  wind  profile  estimates  .-Xlthough  computations  can  be.  and  ai  e.  carried  out  for  the  nonh 
component  of  the  thermal  wind  profile,  the  magniiiKles  are  i_\pic;illy  uiihin  the  ‘‘noise"  of  the 
asailable  data  and  should  not  be  considered  very  meaningful. 

In  order  to  reflect  an  operational  situation  with  sunie  realism,  the  "exact"  temperature 
gradient  profiles  in  Fig.  2.2-.'  were  not  used  in  most  of  the  following  results.  Instead,  reasonable 
piecewise  linear  approximations  to  the  two  compi'iients  were  used  which  are  representative  of 
what  might  be  generated  from,  for  example,  forecasts  of  the  temperature  fields  in  the  region. 
Figure  2.2-4  shows  three  different  tipproximatii'ns  ui  the  (larger)  north  component  and  the  sin¬ 
gle  approximation  used  for  the  (smaller)  east  component,  each  superimposed  over  the  corre¬ 
sponding  data-derived  profile.  Sensitivity  of  the  algorithm  to  the  assumed  profile  will  be 
indicated  by  results  based  on  these  three  north  component  approximtiiions  presented  in  the 
next  subsection. 

2.2.3  Baseline  Optimal  Estimation  Results 

The  primary  purpose  of  this  section  is  to  establish  the  baseline  case  on  which  the  sensi¬ 
tivity  results  presented  in  the  next  subsection  arc  based.  Of  almost  equal  importance  are  the 
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Figure  2.2-4  .Approximate  1  cmperature  Gradient  Profiles 


resLilts  conceining  dependence  ot  the  ne\\  opiimal  estinuilu'ii  ulenmlim  on  three  tundamental 
selections 

•  Choice  of  the  grid  from  which  the  ladiance  gradient  data  are  derived 

•  Choice  of  the  HIRS  2  channels  (i.e..  wavelengths)  used  by  the  data  processor 

•  Choice  of  the  initial  guess  of  the  temperature  gradient  profile. 

Our  baseline  case  will  he  basei.1  on  the  following  selections; 

•  Radiance  gradients  computei.1  from  the  7  \  7  grids  (as  discussed  in  Sec¬ 
tion  2.1) 

•  Channels  3,  4.  5.  o,  15.  and  lb 

•  The  interpolated  temperature  gradient  profiles  in  lag  2.2-3. 

Selection  of  the  "  grid  ividiance  gradient  has  been  (.li'Ciisscd  in  Section  2.1.  The  si.\  IR  chan¬ 
nels  were  chosen  because  their  kwels  of  peak  energ\  contributions  a|)pear  within  the 
troposphere  and  are  less  likel>  to  lie  influenced  b>  earth  surface  effects.  Finally,  while  using  the 
interpolated  temperature  gradient  profile  may  not  be  realistic  operationallv.  its  use  here  as  the 
baseline  permits  more  accurate  evaluation  of  the  algorithm's  sensiti\ity  with  respect  to  model 
parameters.  Detailed  investigations  of  the  influence  of  the  choice  of  appro.ximate  temperature 
gradient  profiles  and  of  channels  are  summarized  in  this  subsection. 

Results  for  The  Baseline  Case  are  summarized  in  Table  2.2-2.  In  all  of  the  tables  in  this 
'ubsection.  estimates  entered  on  a  line  correspon(.ling  to  a  specific  pressure  are  of  the  difference 
between  the  geostrophic  wind  component  at  that  level  and  that  at  the  reference  level  (850  mb  in 
all  cases).  This  table  displays  both  the  initial  theinitil  wind  estimate  which  is  consistent  with  the 


Table  2.2-2  Results  of  the  Baseline  Case"' 
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initial  teniperatLire  gradi-iit  profile  and  the  final  estimtite  ^vliieit  re-'iilts  from  processing  the  sat¬ 
ellite  dtita.  Note  the  significant  difference  in  the  relative  magnitnde  of  the  two  wind  components 
(in  both  tile  initial  and  final  estimates);  this  is  consistent  with  the  etirlier  comments  regarding 
the  nuich  smaller  temperature  deriviiti'.e  in  the  east-west  direction.  Note  also  that  there  is  sig¬ 
nificant  difference  between  the  initial  tiiKl  final  tliermal  w  ind  estimates.  This  is  the  first  clear  in- 
dicahon  that  the  satellite  radiance  data  are.  in  some  sense,  inconsistent  with  the  initia. 
temperature  analysis  fiehl. 

Dependence  on  Temperature  (iradient  Profile 

.A.S  noted  etirlier,  an  operational  algorithm  cannot  expect  to  ha\e  a  complete  tempera¬ 
ture  gradient  profile  derived  from  iiulepeiulent  soundings  a^allable  at  forecasting  times.  In  or¬ 
der  to  study  the  effect  of  using  an  ap(iro.\imate  profile,  the  four  piecew ise-linear 
approximations  to  the  interpolated  profile  which  were  shown  in  Fig.  2.2-4  are  used  in  the  algo¬ 
rithm  (T.\  in  the  figure  represents  an  approximation  to  the  east  derivative  and  T'l  I.  T\’2  and 
T't  s  are  approxinnitioiis  to  the  north  derivative.)  l  ables  2.2-3  and  2.2-4  respectively  contain  the 
initial  and  final  thermtil  wind  estimates  based  vin  etich  of  these  approximations  (as  well  as  the 
baseline  results  for  comparison), 

.■\s  would  be  ex[iectevl,  the  goovi  piecewise  linear  approximations  to  the  interpolated 
metisurements  (T.\  and  iV3)  vield  results  which  are  close  to  those  obtained  using  the  exact 
measured  profile  at  most  levels  [he  sign  differences  which  tippear  in  the  nonh-component 
estimates  do  not  represent  significxint  vleviaiiv>ns  hcciiiise  the  mtigmtudes  are  very  small.  The 
other  two  east-comiioneiit  approximtitions  result  in  estimated  profile  ,  '  represented  by 

Table  2.2-3  Initial  Tliermal  Wind  Estimates''* 


<  I  I  W'''ikI  spectls  iii  m  s 
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Table  2.2-4  Final  Thermal  Wind  Kstiinates'' 
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^allle^  ac  distinct  lc\c.  U  i  winch  are  M)niewhat  ftirihcr  trom  ilie  Ba'-eline  esliinate.  The  iinpor- 
uinr  point  is  that  ulthoii'Ji  the  inituii  profiles  ( orrespondm-^  to  these  opiiriinniations  arc  substantially 
different  than  the  Baseline  initiai  profile,  the  final  esinnate^  are  minJi  ( loser  to  those  Jor  the  Baseline. 
That  is.  there  is  a  coinerccnce  tovsard  a  coninion  protile  tis  a  result  of  processing  the  satellite  ra¬ 
diance  gradient  data.  I, See  alsi>  the  Hxecufive  .Sun)n);)r>  and  f'ig.  BS-2). 

Dependence  on  Choice  of  Channels 

■As  indicated  ahose.  the  si.\  IR  channels  in  l  ahle  2  1-1  were  chosen  because  of  the  lo¬ 
cation  of  their  peak  absorption  bands  and  their  irelaiive)  lack  of  sensitivity  to  surface  radiance. 
The  purpose  of  this  brief  section  is  to  present  one  set  of  results  based  on  subsets  of  those  si.x 
channels.  Again,  only  the  east-components  of  the  thermal  wind  estimates  are  considered. 
.Appro.ximation  TS  2  was  used  for  the  temperature  derivati\e  profile  in  this  study  to  reflect  a 
non-ideal  operating  condition.  Talde  2  2-.'s  contains  the  ti\e-level  profiles  for  six  selections  of 
channel  sets  along  with  the  si.x-chitnnel  results  from  lable  2.2-4 

Comparing  Case  11  with  the  all-channels  result  indicates  that  the  omitted  channel  4  data 
can  be  eliminated  with  little  effect  /or  this  case.  On  the  other  hand,  while  dropping  channel  6 
(Case  1)  produces  a  large  change  in  the  estimated  profile,  dropping  both  channels  6  and  15 
(Case  III)  again  produces  a  result  close  to  the  all-channels  result.  Case  IV,  another  set  of  four 
channels,  is  significantly  poorer  than  Case  111.  Cases  V  and  VI,  which  use  only  three  channels, 
produce  equally  poor  estimates. 
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Tal)k“  2.2-5  Tliennal  Hast  NNiiul  Profile  Kstimate 
Dependence  on  Channel  Selection'" 
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In  e.xamining  thcNc  tcMilt'  it  >hiHikl  Ik-  i'cnic.‘inhci'cd  that  the  iintlei'l>  ing  idea  is  to  find  a 
linear  combination  ol  tiinctittns  ithc  kernel  lunctlon^.  K. )  •'Och  that  the  statistical  optimization 
criterion  is  minimized.  7  hereloi  e.  using  feAver  ol  tlt^tse  Itinctions  can  onl\  increase  the  e.xpected 
square  error.  From  the  point  ot  \iess  of  looking  for  an  operational  implementation  ot  the  algo¬ 
rithm.  using  four  chiinnels  (iiistetid  of  six)  otfers  an  insignificant  saving  in  contputer  time  and 
dtita  handling  requirements,  so  is  not  ;i  critical  design  issue.  It  would  be  possible  to  investigate 
the  contributions  of  each  chtinnel  over  a  stiitistictilly  significant  data  set  and  make 
recommendations  regarding  the  best  choice  for  channel  wtivclengths.  but  such  an  effon  is  be¬ 
yond  the  scope  of  this  research  contivict. 

2.3  OHRING  S  METHOD 

The  approach  to  wind  shettr  profile  estimation  proposed  by  Ohring.  Neeman.  and  Dun- 
ctin  (Ref.  1)  was  discussed  in  Ref.  2  as  a  point  of  departure  for  the  new  optimal  estimation 
approach  under  investigation.  In  order  to  compare  results  of  the  two  approaches,  the  basic  algo¬ 
rithm  in  Ref.  1  (hereafter  refeiTed  to  as  Ohring's  method)  is  outlined  in  this  section.  These 
equations  will  be  used  again  later  while  discussing  sensitivity  of  the  two  algorithms. 

The  estimate  of  the  wind  shear  (east  component)  between  levels  z,  and  z,»i  is  based  on  a 
weighted  sum  of  metisured  (actually,  estimtited)  radiance  gradients  just  as  is  the  optimal  estima¬ 
tor  developed  in  this  research, 

U|(i)  =  -(/?/f)  c,‘  R, 


(2.3-1) 


where  the  sub>cript  L  indicates  an  estinitiie  of  wind  shcai'  mirr  one  specific  atmospheric  layer, 
and 


Rs  = 


(2.3-2j 


(m  =  number  of  radiance  measurement  cliannels).  and  tlie  i  subscript  in  c,  emphasizes  that  the 
optimal  weights  depend  on  tlie  layer  through  which  the  tliermal  windsliear  is  estimated.  The 
vector  of  weialits.  c, .  is  calcuhited  to  minimize  the  following  lunction: 


(c  )  = 


K,(C) 


c) 


dC 


(2.3-3) 


where 


w,i2) 


1  •  C  ■  /|.i 

0  oiiieiwise 


(2.3-4) 


The  Kj,  j  =  1 . m.  are  the  same  kernel  functions  used  in  the  new  optimtil  estimation  approach. 

This  corresponds  to  "Method  1"  in  Ref.  1.  .-\lthough  four  other  minimiz.ation  criteria  were 
evaluated  in  that  reference,  there  was  not  a  sufficient  improvement  using  the  more  complex  cri¬ 
teria  which  warranted  their  intplementation  for  this  comparison.  The  optimal  weighting  vector 
is  calculated  using. 


where 


and 


c,  -  M  '  b, 


i: 


M  =  K(c)  KfO'dc 


•  ^  I 


k(c)  dC.  i  =  I . n 


(2.3-5) 

(2.3-6) 

(2.3-7) 


Two  results  using  Ohring's  method  are  presented,  one  using  all  six  HlRS/2  channels,  and 
the  second  using  a  subset  containing  channels  3,  4,  and  15.  While  these  tv\o  cases  will  be  com¬ 
pared  with  results  of  the  Optimal  Estimation  Approach  presented  in  the  previous  subsection, 
they  certainly  do  not  constitute  a  complete  evaluation  of  Ohring’s  method  for  the  sample  data 
set.  As  will  be  shown,  kernel  functions  for  the  six  HIRS/2  channels  are  far  from  ideal  for 
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appioxiniating  the  pLll^e•^l1aped  ucigliting  tunctions  (Fa].  2.3-4).  This  condition  leads  to  severe 
numerical  problems  uhen  implementing  Ohring’s  method  in  a  straightforward  manner. 

Using  the  radiance  grtidieiits  computed  from  ^  ,\  7  grids  (ax  discussed  in  Section  2.1) 
and  all  si.x  HIRS  2  channeK  resulted  in  the  extremely  poor  estimates  of  the  East  thermal  wind 
within  each  layer  ^hown  in  1  able  2.3- 1 .  Fhe  reason  for  this  poor  result  is  that  the  least-squares 
problem  (see  Fa].  2.3-3)  leads  to  a  minimi/ing  solution  vector,  cx ,  for  layer  i  which  has  a  very 
laree  norm  becau  .e  o;  .he  piK)rly  conditioned  M  matrix  (FT).  2.3-0).  This  poor  conditioning  is  a 
direct  result  of  the  dcpcndcix  e  timong  the  kernel  functions  (see  -Vppendix  C).  That  is,  the  kernel 
ftinctions  are  too  simikir  so  th;it  there  exist  linear  combinations  of  them  which  are  nearly  identi¬ 
cally  zero  on  each  interval  of  interest 

That  the  least-squiires  [iroblem  has  been  solved  ciinectiv  is  demonstrated  in  Figs.  2.3-1 
throuiih  2.3-5  which  show  the  (ideal)  puKc  function  and  the  approximating  curve  defined  by  the 
"optimal "  linear  combination  of  the  six  kernel  functunis  (solid  cuiAes)  for  each  of  the  five  lay¬ 
ers  within  which  the  thermal  wind  is  estimated  Clearly,  the  approximating  functions  have  large 
non-zero  portions  outside  of  the  estimated  la>er.  a  direct  coiiseciuence  of  the  enforced  minimi¬ 
zation  of  the  criterion  function. 

For  the  three-clutnnel  ctise.  the  iippro.ximating  functions  (ilashed  lines  in  Figs.  2.3-1  — 
2.3-5)  while  not  as  "shaip  '  in  the  region  of  the  ideal  pulses,  avoid  much  of  the  large  side-lobe 
oscillations.  (See,  tor  example,  the  much  smoother  curves  for  the  two  layers  above  300  mb.) 
But.  as  shown  in  the  right  hand  column  of  f  able  2.3-1,  this  three-channel  case  results  in  more 
reasonable  estimates  of  the  Fast  thermal  wind  in  each  layer. 


Table  2.3-1  Ohring’s  Method  Results:  East  Component 
of  Thermal  Wind 


L.A'i  FR 

0  CFf.-XNNELS 

3  CU.-\NNELS 

,S50  -  ^’OO  mb 

-5.0 

3.0 

7()0  -  s()()  mb 

-132.1 

8.1 

500  —  300  mb 

172.1 

13.4 

300  -  200  mb 

109.1 

h.'-) 

200  -  100  mb 

-431.4 

5.4 
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Figure  2.3-3  Weigtning  FLiiK'tK)n  Approxiniations  on  500  mb  to  300  mb 
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Figure  2.3-4  Weighting  Function  Appro.ximations  on  300  mb  to  200  mb 
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Figure  2.3-5  Weighting  Function  Appro.\iniations  on  200  mb  to  100  mb 


Table  2.3-2  contains  tlie  calculated  coefficients  for  tuo  typictil  layers  for  the  two  cases 
and  Table  2  3-3  shows  the  East  component  estimate  contributions  resulting  from  these  coeffi¬ 
cients.  Clearly,  the  large  coefficients  in  the  si.x-channel  case  lead  to  sums  of  positive  and 
negative  terms  whose  magnitudes  are  much  larger  than  the  eNtimated  quantities.  This  effect  is 
especially  clear  in  channels  4.  5  and  6,  a  direct  result  of  the  similar  kernel  functions  for  these 
channels  (see  Appendi.x  C). 

Table  2.3-2  Optiinul  Coefficients  from  Ohring’s  Method  (East  Component) 


CHANNEL 

LA\ER  1 

!S5()  mb 
6  channels 

-  700  mb 

3  channels 

200  mb  - 
6  channels 

-  100  mb 

3  channels 

3 

-557 

-213 

4951 

-1381 

4 

3963 

301 

-4,05x10-’ 

193 

5 

-9730 

— 

6.56x10-’ 

6 

4986 

-2.04x10’ 

— 

15 

-5.50x10-’ 

-  1.85x10-’ 

-4.94x1 0  = 

-1987 

16 

4.37x105 

-8.19x105 
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Table  2.3-3  Optimal  Ohrin^’s  Metliotl  Contributions  (East  Component) 


CH,‘\NNEL 

L.A\ER  1 

X50  mb  - 

6  channels 

-  700  mb 

3  channels 

200  mb  — 
0  channels 

100  mb 

3  channels 

J) 

3.4 

1.3 

-20.7 

8.3 

4 

-7(1.3 

-5.8 

780.9 

-3.7 

5 

451.5 

--- 

-3044.8 

— 

6 

-40(1.2 

— 

l(i(i2.2 

— 

15 

22.4 

--- 

200.7 

— 

Ih 

0.4 

7  5 

-0.7 

0.8 

TOT.AL 

-.TO 

3  0 

-431.4 

5.4 

Finally,  Table  2.3-4  lists  the  estimates  tor  both  components  ot  the  thermal  wind  profile 
produced  by  Ohring's  method  (using  only  chantiels  3.  4  and  15)  along  with  those  generated  by 
the  Optimal  Estimation  t  echnique.  (The  East  component  of  this  profile  was  presented  earlier  in 
Table  2.2-5.  case  \  )  Note  that  the  values  in  Table  2.3-4  are  the  estimated  thermal  wind 
between  the  indicated  level  and  the  reference.  ■''50  mb.  level;  the  east  component  values,  for  e,\- 
ample,  are  obtained  from  the  results  in  fable  2.3- 1  by  summing  the  layer  estimates  up  to  the  in¬ 
dicated  pressure  level.  It  should  be  emphasized  that  tiiere  is  no  fundamental  difference  between 


Table  2.3-4  Coinpurisoii  of  Optimal  Estimation  and  Ohring’s 
Method  Estimates  (Channels  3.  4.  and  15) 


WIND  PROFILE  (m/ 

s) 

LEVEL 

East 

North 

Ohring 

Opt.  Est. 

Ohring 

Opt.  Esi. 

700  mb 

3.0 

4.4 

-5.2 

-0.4 

500  mb 

1  1.0 

10.3 

1.0 

-0.6 

300  mb 

24.5 

15.3 

21J 

-0.0 

200  mb 

31.4 

15.9 

39.0 

0.4 

100  mb 

36.8 

9.9 

33.0 

0.5 
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the  two  formulations  (i.e..  layer-by-layer  estimate'.  \s.  summed  layers);  Ohring's  method  and 
the  Optimal  Estimation  approach  have  much  different  characteristics,  but  not  because  of  this 
feature  of  their  fomuilation. 

To  summarize,  although  Ohring's  method  is  susceptible  to  numerical  difficulties  due  to 
dependent  kernel  functions,  it  could  be  stabilized  by  using  a  singular  value  decomposition 
procedure  or  any  similarly  robust  algorithm  for  ill-conditioned  least-squares  problems  (Ref.  3). 
The  impotiant  difference  between  Ohring’s  method  and  the  Optimal  Estimation  approach  is 
that  because  the  former  method  is  based  on  obtaining  a  least-squares  fit  to  the  pulse  functions 
(Eq.  2.3-4),  its  estimated  profiles  can  never  achieve  a  smaller  mean-square  estimation  error 
than  those  of  the  latter  approach  over  a  representative  sample  of  cases. 


2.4  SENSITIVm'  EVALL.UTONS 

The  objective  of  this  sensitivity  evaluation  is  to  demonstrate  the  dependence  of  thermal 
wind  profile  estimates  on  algorithm  data  and  parameters.  The  first  subsection  deals  only  with 
the  new  optimal  estimation  algorithm  and  uses  direct  penurbation  of  the  data  and  parameters  to 
evaluate  sensitivity.  In  the  second  subsection,  analytical  e.xpressions  derived  in  Appendi.x  F  for 
logarithmic  derivatives  of  the  estimates  of  both  the  optimal  estimation  and  Ohring  algorithms 
will  be  used.  These  local  results  describe  sensitivity  with  respect  to  small  changes  in  the  kernel 
functions  and  provide  a  basis  for  two  interesting  comparisons:  fii  st.  with  the  large-scale  penur¬ 
bation  results  for  kernel  functions  presented  in  the  first  subsection,  and  second,  between  the  lo¬ 
cal  sensitivities  of  the  optimal  estimation  and  Ohring  algorithms. 

Results  in  this  section  augment  the  investigations  reported  in  Section  2.2.2  into  depen¬ 
dence  of  the  optimal  estimation  algorithm  on  the  assumed  temperature  gradient  profile  and  on 
the  choice  of  HlRS/2  channels.  Emphasis  in  this  section  is  on  perturbation  of  scalar  parameters 
instead  of  those  more  comple.x  algorithm  variations,  I  he  kernel  function  sensitivities  reported 
here  are.  in  fact,  also  based  on  functional  perturbations  although  those  perturbations  are  de¬ 
fined  by  scalar  parameters. 
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2.4.1  Perturbations  of  Data  and  Parameters 


All  penurbaticns  prc'-.eracd  in  this  section  are  depanures  from  the  data  and  parameters 
Used  in  the  baseline  ease  described  in  Section  2.2.2  for  the  optimal  estimation  algorithm.  Sensi¬ 
tivities  with  respect  to  the  following  data  and  parameters  have  been  investigated; 

•  Temperature  gradient  profile 

•  Radiance  gradients 

•  Radiance  gradient  error  covariance  matri.x 

•  Kernel  functions. 

Sensitivities  were  evaluated  separately  for  each  item,  then  a  single  run  was  made  with  all  of 
these  items  penurbed  to  assess  the  combined  effects  in  at  least  one  case.  Each  of  the 
investigations  required  a  slightly  different  perturbation  technique;  these  are  discussed  individu¬ 
ally  along  with  each  result  in  the  remainder  of  this  subsection. 

Results  are  e.xpressed  in  normalized  form;  fractional  (or  percentage)  change  in  the  esti¬ 
mate  is  plotted  against  fractional  change  in  the  parameter.  One  can  define  a  sca'ar  sensitivity, 
based  on  a  linear  model  by 

S,,  =  |Au,/u;’|  /  lApj/p'j’l  (2.4-1) 

where  the  superscript "  represents  baseline  values.  Values  of  S,,  will  be  used  in  the  following  dis¬ 
cussion  when  comparing  senvitivities  with  respect  to  different  parameters.  Finally,  results  are 
presented  only  for  the  effect  of  penurbations  on  the  zonal  (i.e..  east)  component  of  the  thennal 
wind  estimate.  This  choice  was  made  because  the  meridional  (north)  component  of  the  thennal 
wind  estimate  was  too  small  to  be  significant. 

i)  Temperature  Gradient  ProFde  Errors 

Zonal  wind  estimates  depend  on  the  meridional  component  of  the  temperature  gradient 
profile  (Eqs.  2.2-5  and  6).  Penurbations  were  first  applied  to  the  temperature  gradients 
developed  from  the  NMC  analysis  fields  (see  Section  2.2.2)  at  two  individual  levels  (850  mb 
and  500  mb)  and  then  simultaneously  at  all  si.x  levels.  For  each  case,  a  continuous  profile  was 
generated  by  linear  interpolation  between  the  six  levels,  just  as  was  done  to  develop  the  baseline 
profile  in  Fig.  2.2-3.  Because  of  this  procedure,  perturbations  at  one  level  cause  profile  pertur¬ 
bations  only  between  adjacent  levels;  from  1013  mb  (surface)  to  700  mb  for  the  850  mb  pertur¬ 
bation  and  from  700  mb  to  300  mb  for  the  500  mb  perturbation. 
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Two  types  of  peniirbations  were  used:  additive  and  nuiltiplicative.  Additive  penurba- 
tions  are  of  the  form 


T^(/,)  =  +  r<|  r;(/,)l  (2.4-2) 

where  represents  the  panial  derivative  of  temperature  with  respect  to  y  (the  north  coordi¬ 
nate).  The  perturbation  temi  in  Eq.  2.4-2  is  ahvays positive  and  equal  to  the  selected  fraction,  a  , 
of  T^(Zi).  Additive  penurbations  are  characteristic  of  bias  or  "drift"  instrumentation  errors  or 
data  processing  algorithm  effects.  Multiplicative  perturbations  are  given  by 

T^(7.,)  =  (I  -I-  a)  Ty(/,)  (2.4-3) 


where  a  is  the  selected  fractional  peiturbation.  Multiplicati\e  perturbations  are  characteristic 
of  scale-factor  instrumentation  errors  or  data  processing  algorithm  effects.  That  these  two 
forms  lead  to  significantly  different  temperature  profiles  is  illustrated  in  Fig.  2.4-1  which  shows 
the  nominal  and  both  perturbations  of  the  meridional  "‘mperature  gradient  profile  based  on  a 
40'^c  perturbation  at  all  levels. 


f  -  !  4  -  -s 


Figure  2.4-1  Nominal  and  Perturbed  Meridional  Temperature 
Gradient  Profiles 
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Results  toi'  additi\e  [■'ei  tLirhations  of  l()''c.  2()‘^f  and  in  the  entire  meridional  temper¬ 
ature  dernative  profile  on  iliermal  wind  estimates  within  three  layers  are  summarized  in 
Table  2  4-1  The  radical  i^litference  between  sensitivities  of  estimates  in  layers  below  200  mb 
and  that  in  the  200  mb  to  100  mb  Itiver  can  be  understood  through  careful  examination  of  the 
estimation  equations  [  ig  2  2-3  re\eals  that  the  nominal  temperature  gradient  profile  reverses 
sign  \Mthin  the  highK  sensiti\e  layer,  and  that  sign  reversal  leads  to  a  penurbation  term  which 
does  not  exist  in  other  layers.  Within  the  other  layers  sensitivity  actually  decreases  with  increas¬ 
ing  penurbation.  almost  sLirely  due  to  numerical  effects. 

•Additive  penuibations  of  the  meridional  temperature  gradient  at  only  one  le\el  pro¬ 
duced  the  results  shown  in  Fig.  2  4-2.  (perturbing  the  gradient  at  850  mb)  and  Fig.  2.4-3.  (for 
the  gradient  at  500  mb).  The  positise  perturbations  decreased  the  themial  wind  estimates  in 
adjacent  layers,  but  increased  estimates  within  other  layers  in  every  case.  Values  of  the  sensitiv¬ 
ity  coefficient.  S.  range  between  0  20  and  0.56  for  these  perturbations.  That  is.  the  fractional 
cliange  in  the  estimates  were  smaller  than  the  fractional  changes  in  the  data  by  44‘^f  to  SO'^T, 

Multiplicative  perturbations  produced  different  sensitivity  characteristics,  but  similar 
sensitivity  coefficient  salues.  F’erturbing  the  temperature  gradient  values  at  all  levels  has  almost 
no  effect  on  the  estimates.  This  is  a  direct  consequence  of  the  approximate  linearity  of  the  esti¬ 
mates  in  the  temperature  gradient  profile  for  the  radiance  gradient  sensitivity  model  used 
during  this  demonstration.  The  linearity  is  due  to  the  fact  that,  in  Eq.  2.2-4.  the  diagonal  ele¬ 
ments  of  the  second  matrix  (dependent  onf^)  dominate  the  gamma  matrix  diagonals. 

When  multiplicative  pentirbations  are  applied  to  \alues  at  only  selected  levels,  larger, 
but  still  rather  small  sensitivity  is  evidenced.  Figure  2.4-4  summarizes  the  results  for  perturba¬ 
tion-,  111  the  meridional  temperature  derivative  at  ,850  mb.  Tlie  S  values  for  this  case.  0.1  and 
0.39,  are  typical  of  those  for  multiplicative  penuibations  at  other  levels. 

Table  2.4-1  Fractional  Changes  in  Layer  Estimates 
Due  to  Additive  Perturbations 
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Figure  2.4-2  Sensitivity  of  Thermal  Wind  Kstimates  to  a  850  mb  Temperature 
Gradient  Additive  Perturbation 


ii)  Radiance  Gradient  Perturbations 

The  optimal  themial  wind  estimation  algorithm  is  a  linear  mapping  of  the  measured  ra¬ 
diance  gradients  into  thermal  wind  estimates,  rherctoic.  perturbation  of  all  radiance  gradient 
channels  according  to 


VR  =  (VR)„(1  +  a)  (2.4-4) 

for  any  a  yields  a  thermal  wind  estimate  scaled  by  the  same  ( 1  +  «  )  factor.  (In  effect,  S  =  1  for 
this  case.) 
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Figure  2.4-3  Sensitivity  of  Thermal  Wind  Estimates  to  a  500  mb  Temperature 
Gradient  Additive  Penurbation 


More  interesting  results  are  obtained  by  considering  penurbaiions  in  each  channel  sepa¬ 
rately.  In  fact,  it  is  easy  to  show  that,  for  the  thermal  wind  estimate  within  any  layer  resulting 
from  perturbing  the  j'^  radiance  gradient, 

^  =  n  (a,/u")  (2.4-5) 

u 

where  aj  is  the  contribution  to  the  unperturbed  estimate,  u”,  by  )^Rj/3y.  Table  2.4-2  lists  the  fac¬ 
tor  Sj  =  aj/u"  in  Eq  2.4-5  for  each  of  the  si,\  channels  used  in  the  demonstration.  Channels  5, 
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Figure  2.4-4  Sensitivity  of  Thermal  Wind  Estimates  to  a  S5()  mb  Temperature 
Gradient  Multiplicative  Perturbation 

6  and  15  are  clearly  the  major  contributors  to  the  estimates,  therefore  exhibiting  the  highest 
sensitivity. 

iii)  Radiance  Gradient  Error  Variance  Perturbations 

These  parameters,  the  diagonal  elements  of  the  radiance  gradient  estimation  error 
covariance  matrix,  F  in  Eq.  2.2-4.  were  included  in  the  perturbation  analysis  because  their 
“correct”  values  would  vary  considerably  in  an  operational  system  due  to  variations  in  the 
radiance-measuring  instrument,  the  geometry  of  radiance  measurements,  and  possibly  atmo¬ 
spheric  conditions.  If  these  variations  were  modeled  accurately,  significant  changes  in  the  error 
variances  could  result.  (Of  course,  an  implementation  of  the  algorithm  could  ignore  the 
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Table  2.4-2  Sensitivities  with  Respect  to  Radiance 
Gradient  Components 


CHANNEL 

S 

3 

0.0029 

4 

0.0216 

5 

0.2948 

6 

0.5220 

15 

0.1645 

16 

0.0000 

variations  and  use  a  conservative,  but  incoriect.  large  variance  in  all  conditions.  Consideration 
of  this  alternative  is  beyond  the  scope  of  the  current  ^tudy.) 

For  simplicity,  all  diagonal  components  were  penurbed  by  an  identical  scale  factor.  This 
i^  representative  of  the  effects  of  altered  radiance  data  location  geometiy,  but  with  the  individu¬ 
al  channel  error  variances  either  remaining  fi.xed  or  scaling  equally.  Results  for  this  investiga¬ 
tion  are  presented  in  Fig.  2.4-5.  The  relatively  small  setisitivity  of  estimates  to  these  parameters 
is  dramatically  indicated  by  the  logarithmic  scale  on  the  abscissa.  Note  that  for  scaling  by  factors 
less  than  1.0.  the  sensitivity  remains  zero.  The  e.xtremely  small  sensitivity  for  factors  less  than 
100  is  a  consequence  of  the  P  matrix  (F.q.  2.2-4)  being  dominated  by  the  term  rather  than 
the  r  matrix. 


iv)  Kernel  Function  Perturbations 

Kernel  function  penurbations  are  the  most  critical  of  those  considered  because  of  the 
known  sensitivity  of  the  kemei  functions  themselves  to  variations  in  the  temperature  profile  on 
which  they  are  based.  Because  an  operational  implementation  would,  of  necessity,  use  at  best  a 
forecast  of  the  temperature  profile  at  the  time  of  data  collection,  this  known  sensitivity  could 
lead  to  significant  thermal  wind  estimation  error. 

Kernel  function  sensitivity  to  temperature  variation  was  in\estigated  analytically  by  dif¬ 
ferentiation  and  numerically  by  making  finite  penurbations  in  the  temperature  data.  From  the 
kernel  function  formulas  found  in  Appendix  G,  one  can  show  that  (approximately): 


Sij(z)  = 


T(7) 

K,(z) 


fiKj 

la 


(z) 


(2.4-6) 
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Figure  2.4-5  Sensitivity  of  Thermal  Wind  Estimates  to  Radiance 
Gradient  Error  Variance 


where 


C2  =  hc/k  =  1.439  cm  deg 


(2.4-7) 


(h  and  k  are  the  Planck  and  Boltzman  constants  and  c  is  the  speed  of  light),  and  T(z)  is  the  ac 
tual  temperature  at  the  level,  p(z),  of  interest  and  is  the  wave  number  (in  cm’’)  of  channel  j. 
The  maximum  error  in  Eq.  2,4-6  is  under  6%  and  occurs  at  that  z  for  which  the  temperature  is  a 
minimum  and  for  the  minimum  wave  number  channel  (j  =  3).  Applying  Eq.  2.4-6  yields  the 
(infinitesimal)  sensitivities  summarized  in  Table  2.4-3.  Recall  that  Stj  =  2.0,  for  example, 
means  that  a  1%  change  in  T(z)  would  cause  a  2%  change  in  Kj(z).  Thus,  the  sensitivities  in 
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Table  2.4-3  Sensitivity  of  Kernel  Functions  to  Temperature 


CHANNELS 

SENSITIVITY  (S,j)  RANGE 

3.  4.  5.  6 

1.5  <  Stj  <  3.0 

(small  wavenumbers) 

15.  16 

9.0  <  St,  <  13.5 

(large  wavenumbers) 

Table  2.4-3  are  indeed  large.  However,  a  I'^c  change  in  temperature  is.  in  fact,  quite  large  rela¬ 
tive  to  typical  temperature  uncenainties  for  nominal  temperatures  between  210"  K  and  289"  K. 
If  a  profile  were  in  error  by  1"  K  over  its  entire  length  (i.e..  less  than  ().5‘^c),  the  percentage 
change  in  Kj(z)  would  be  on  the  order  of  O.S'^r  and  for  the  low  and  high  wave-number  chan¬ 
nels,  respectively. 


For  the  numerical  investigation  of  kernel  function  sensitivity  to  temperature,  equal  (posi¬ 
tive)  percentage  changes  were  made  to  the  entire  temperature  profile,  the  resulting  perturbed 
kernel  functions  computed,  and  the  average  percentage  deviation  of  Ki(z)  from  its  nominal  val¬ 
ues  between  z  =  0  (the  reference  pressure  level,  p  =  1013  mb)  and  the  top  of  the  atmosphere  was 
calculated.  Figure  2.4-6  contains  plots  of  this  average  sensitivity  for  channels  typical  of  the 
smaller  wave  numbers  (channel  3)  and  larger  wave  numbers  (channel  15).  The  dramatic  differ¬ 
ence  in  the  behavior  of  these  two  cases  is  due  to  the  effects  of  the  e.xponential.  exp  C2»'j/T(z)|. 
which  appears  in  the  Planck  function  factor  of  the  kernel  function  formulas  (see  Appendi.x  F). 
Clearly,  the  extreme  sensitivity  of  the  h.gher  wave  number  channels  already  apparent  in  the  in¬ 
finitesimal  sensitivities  in  Table  2.4-3  is  even  more  evident  in  the  finite  perturbation  results.  As 
a  result,  an  operational  implementation  of  the  optimal  estimation  algorithm  examined  in  this 
project  should  attempt  to  use  relatively  smaller  wave  number  (i.e.,  longer  wave  length) 
channels. 


Finally,  the  impact  of  unifonn  percentage  changes  in  the  kernel  functions  on  the  thennal 
wind  estimates  is  shown  in  Fig.  2.4-7.  Estimates  are  most  sensitive  to  channel  15  for  which  Sj  = 
0.5  for  up  to  a  26%  change  in  the  kernel  function  (at  all  values  of  z).  Perturbations  in  all  chan¬ 
nels,  which  would  be  the  case  when  perturbations  were  due  to  a  temperature  profile  error,  lead 
to  St  =  0.8. 
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Figure  2.4-6  Average  Sensitivity  of  the  Kernel  Functions  to  the 
Temperature  Profile 


The  investigation  presented  here  did  not  pursue  the  obvious  next  step,  calculating  esti¬ 
mate  sensitivity  due  directly  to  temperature  profile  peiturbation,  because  realistic  perturbation 
profiles  were  not  available  during  conduct  of  this  study.  Because  of  the  relatively  large  sensitiv¬ 
ity  of  the  higher  wave  number  channels,  it  is  extremely  important  to  assess  this  sensitivity.  Real¬ 
istic  temperature  profile  perturbations,  rather  than  arbitraiy  constant  percentage  perturbations, 
should  be  used  in  the  study  to  avoid  raising  unrealistic  concerns  due  to  improbable  or  impossi¬ 
ble  parameter  perturbations. 
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Figure  2.4-7  Sensitivity  of  Thermal  Wind  Estimates  to  the 
Kernel  Functions 
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v)  .\11-Paranieter  Perturbations 


To  evaluate  the  effect  of  all  paiameters  being  penurbed  together  in  one  simple  way,  ad¬ 
ditive  perturbations  of  the  complete  collection  of  parameters  discussed  in  this  section  were 
made  simultaneously.  Consistent  with  previous  results,  the  percentage  changes  in  thermal  wind 
estimates  across  the  lower  four  layers  were  different  than  that  in  the  topmost  (200  mb  - 
100  mb)  layer.  Table  2.4-4  summarizes  results  for  penurbations  up  to  20'^ c.  Note  that  the  frac¬ 
tional  change  for  the  sum  over  all  layers  is  less  than  that  over  the  first  four.  This  occutred 
because  the  estimate  in  the  top  layer  becomes  more  negative  by  30''T  tlierefore  countering 
much  of  the  positive  (1 1.2'^f)  perturbation  over  the  lower  lasers.  The  sensitivity  (S)  values  are 
included  in  the  table  to  demonstrate  the  non-linear  nature  of  the  penurbation  effects.  That  is,  as 
the  perturbation  magnitude  increased  from  to  the  ratio  of  the  estimate  changes  (over 
all  layers)  to  the  penurbation  magnitude  actuall>  decreased  slightly. 

It  should  be  emphasized  that  although  the  above  simple  e.xperimem  does  provide  some 
useful  insight  into  the  effects  of  multiple  parameter  perturbations,  global  conclusions  regarding 
sensitivity  cannot  be  drawn  from  these  results.  To  obtain  a  more  general  description  of  parame¬ 
ter  sensitivity,  at  least  two  steps  should  be  taken.  First,  the  uniform  perturbation  of  the  kernel 
functions  should  be  replaced  by  a  model  of  typical  temperature  profile  uncenainties.  This 
should  be  a  statistical  model  which  represents  correctly  the  correlation  between  temperature 
penurbations  at  different  levels  in  the  profile.  Second,  a  Monte-C  arlo  procedure  should  be  used 


Table  2.4-4  Fractional  Changes  in  Tlierinai  Wind  *^stiinates 
(All  Parameters  Perturbed) 


PARAMETER  PERTURB, ATION 

LA'tERS 

20^fc 

2-5 

mmm 

Per-Layer 

0,11 

0.24 

Summed 

0.11 

0.24 

■^1 

Sensitivity  (S) 

2.2 

2.4 

2.65 

6 

Per-Layer 

0.30 

0.68 

1.76 

Summed 

0.06 

0.12 

0.20 

Sensitivity  (S) 

1.2 

1.2 

1.0 
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to  sample  penurbations  from  assumed,  reasonable  distributions  of  parameter  uncenainties.  and 
the  distribution  of  the  resulting  estimates  e.xamined.  Parametr'c  models  of  those  distributions 
would  then  provide  a  convenient  way  of  expressing  the  optimal  estimation  algorithm’s  sensitiv¬ 
ity  to  parameter  penurbations  as  a  function  of  assumed  models  of  the  uncenainties  in  those 
parameters. 

2.4.2  Aiialvtic  Sensiiiviiy  Results 


i)  Optimal  Algorithm  Sensitivity 

Equations  F.1-3  -  F.1-9  which  define  the  logarithmic  or  fractional  sensitivity  of  each 
component  of  the  thermal  wind  estimation  vector,  u,  with  respect  to  each  of  the  kernel  func¬ 
tions.  Kj.  were  evaluated  for  the  Baseline  Case  defined  in  Section  2.2.3.  The  results  contained 
one  surprising  but.  on  reflection,  reasonable  result  and  provided  valuable  insight  into  inteipreta- 
tion  of  the  perturbation  results. 

The  surprising  result  is  that  all  rows  of  the  sensitivity  matrix  are  identical,  implying  that 
the  sensitivity  of  the  estimate  is  independent  of  the  (upper)  level,  z.  Table  2.4-5  shows  the  calcu¬ 
lated  sensitivities.  Since  these  numbers  represent  logarithmic  derivatives,  they  should  be  inter¬ 
preted  as  ratios  of  fractional  changes.  For  example. 


Aii, 


=  0.011 


AK3(z) 

M/) 


(2.4-8) 


Table  2.4-5  Kernel  Function  Sensitivity  .Matri.x:  Optimal  Estimate  of 
East  Tliennal  Wind  Component 


LEVELS 

KERNEL  FUNCTION  COMPONENT  (CHANNEL) 

3 

4 

5 

n 

15 

16 

700  mb  ' 
500  mb 
300  mb 
200  mb 
100  mb 

J 

0.01 1 

-0.013 

-0.072 

-0.056 

-0.18 

-0.036 
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where,  on  the  right-hand  side,  the  ratio  is  constant  for  all  values  of  /.  and  the  subscript  indicates 
the  kernel  function  for  channel  3.  Lack  of  dependence  of  the  sensiti\  ity  on  the  level  is  due  to  the 
normalization  in  Eq.  2.4-8  and  is  consistent  with  Lq.  2.2-9  which  sho\\ed  that  the  kernel  func¬ 
tions  affect  the  estimates  only  through  the  scalar,  d.  which  is  independent  of  level. 

The  sign  pattern  (channel  3,  -k.  all  others,  -)  of  the  sensitivities  is  consistent  with  the 
signs  of  the  weights,  c,.  used  to  form  the  optimal  thermal  wind  estimate.  The  relative  magni¬ 
tudes  of  the  temis  are  consistent  with  the  relative  magnitudes  of  the  sensitivities  displayed  in 
Fig.  2.4-7.  The  absolute  magnitudes  of  the  sensitivities  are  much  lower  in  l  able  2.4-5  (by  a  fac¬ 
tor  of  two  to  three).  This  is  due  to  nonlinear  effects  which  drove  the  earlier  perturbation  results 
but  is  not  important  in  the  infinitesimal  case. 

ii)  Ohring  Method  Sensitivity 

The  sensitivity  matri.x  for  Ohring's  method  is  computed  using  Eqs.  F.2-3  -  F.2-9. 
Table  2.4-6  contains  this  sensitivity  m.u.i.x  for  the  three  channel  estimation  case  (see 
Tables  2.3-2,  3.  and  4).  The  contrast  with  the  preceding  case  is  clear.  The  sensitivity  varies  with 
the  upper  limit  of  the  layer  over  which  the  estimate  was  formed  (the  lower  limit  was  850  mb  in 
each  case).  More  importantly,  the  magnitude  of  the  sensitivities  are  much  larger  in  the  (Jhring 
method  sensitivity  matri.x.  This  is  a  significant  result,  provided  it  can  be  shown  to  hold  over  a  va¬ 
riety  of  cases  (i.e.,  data  collection  geometries  and  chtinnel  fretiuencies). 


Table  2.4-6  Kernel  Function  Sensitivity  Matrix:  Ohring  Estimate  of 
East  Tliemial  Wind  Component 


LEVELS 

KERNEL  FUNCTION 
COMPONENT  (CHANNEL) 

3 

4 

15 

700  mb 

-0.35 

n 

-0.33 

500  mb 

-0.16 

WBBM 

-0.16 

300  mb 

-0.087 

0.36 

-0.056 

200  mb 

-0. 1 1 

0.38 

-0.040 

100  mb 

-0.22 

0.54 

-0.059 
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3. 


SUMMARY  AND  RECOMMENDATIONS 


This  Final  Report  has  presented  theoretical  results  and  summarized  a  demonstration  and 
sensitivity  analysis  of  an  optimal  estimation  approach  to  thermal  wind  profile  estimation.  The 
most  important  features  of  ihi.-.  approach  are: 

Direct  estimation  of  the  thermal  winds  from  radiance  data  without  intermediate 
temperature  profile  estimation 

Thermal  wind  estimates  which  are  based  on  error  models  for  all  data  used  by  the 
algorithm 

Thermal  wind  estimates  which  minimize  the  mean  squared  estimation  error. 

The  demonstration  was  based  on  a  single  set  of  radiance  data  and  therefore  did  not  suppon 
assessment  of  the  mean-square  error  properties  of  the  optimal  estimation  algorithm.  Important 
conclusions  which  can  be  drawn  from  the  demonstration  are  described  in  the  following  subsec¬ 
tion.  Recommendations  regarding  fuither  work  are  presented  in  the  concluding  subsection. 


3.1  CONCLUSIONS  FROM  THE  DEMONSTRATION 

Radiance  gradient  calculations  —  Since  the  optimal  estimation  approach  requires  ra¬ 
diance  gradients  for  each  channel  as  input  data,  a  preprocessing  step  is  required  lO  estimate 
these  gradients  from  radiance  measurements  over  a  grid.  During  the  development  and  demon¬ 
stration,  more  time  and  effort  was  spent  than  e.xpecied  on  radiance  gradient  computation  issues 
due  to  the  observed  spatial  variability  of  calculated  radiance  gradients.  Our  conclusion  was  that 
additional  effort  would  be  needed  to  develop  a  robust  algorithm  for  radiance  gradient  calcula¬ 
tions.  A  planer  fit  approach  was  adopted  for  radiance  gradient  estimation  which  can  be  adapt¬ 
able  to  an  operational  environment.  That  adaptation  will  not  be  trivial,  however,  since  it  must 
deal  with  arbitrary  data-array  shapes,  missing  data  points  and  quality  control  to  avoid  cloud-con¬ 
taminated  data. 

Temperature  and  temperature  gradient  proTdes  —  The  sensitivity  results  clearly  indicate 
that  these  are  important  input  data  to  the  algorithm.  For  an  operational  system,  they  would  be 
derived  from  temperature  field  forecasts  or  based  on  persistence  of  a  previous  analysis  field. 
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W'lietlier  or  not  temperature  profile  uncertainties  are  critical  to  successful  application  of  the  al¬ 
gorithm  depends  to  a  large  extent  on  the  form  of  typical  temperature  error  profiles.  As  described 
in  Section  2.4.1(iv),  kernel  functions  are  quite  sensitive  to  bias-type  temperature  profile  pertur¬ 
bations.  That  analysis  should  be  e.xtended  to  realistic  profiles  of  temperature  deviation. 

The  precise  nature  of  the  algorithm  s  dependence  on  the  assumed  temperature  gradient 
profile  was  revealed  by  the  new  result  summarized  in  Eq.  2.2-9.  For  the  estimation  algorithm  im¬ 
plementation  used  during  this  investigation,  the  form  of  the  wind  profile  estimate  is  completely 
detemiined  by  the  initial  temperature  gradient  profile.  Processing  the  satellite  data  will  scale  the 
thermal  wind  corresponding  to  that  profile  so  as  to  best  fit  the  actual  thermal  wind  profile  in  a 
mean-square  error  sense.  Note,  however,  that  this  conclusion  would  change  if  a  vertically  vary¬ 
ing  error  model  for  the  uncertainty  in  the  temperature  gradient  profile  were  included  in  the  im¬ 
plementation  (see  next  paragraph). 

Temperature  gradient  eiror  model  —  It  was  originally  intended  to  develop  and  imple¬ 
ment  a  model  of  the  errors  in  the  assumed  temperature  gradient  profile  for  use  in  the  algorithm. 
By  the  midpoint  of  the  project,  however,  the  extra  time  spent  on  radiance  gradient  calculations 
and  the  desire  for  a  thorough  sensitivity  assessment,  necessitated  delaying  that  development. 
Without  that  error  model,  the  resulting  mean-square  error  estimates  are  unrealistically  low  and 
have  not  been  discussed  in  this  report.  (Typical  values  are  1  m/s  or  less). 

Preliminary  work  was  conducted  on  the  temperature  gradient  error  model;  completing  its 
development  requires  only  selection  of  specific  functional  forms  for  uncenainty  variations  with 
altitude  and  cross-correlations  between  deviations  at  different  levels.  Equations  describing  the 
algorithm  w-ith  this  model  included  have  been  developed.  Their  implementation  is  more  difficult 
than  that  of  the  radiance  gradient  error  model,  but  is  a  >traightfor\vard  task.  Completing  this  im¬ 
plementation  would  be  a  key  part  of  further  work  toward  operational  use  of  the  optimal  estima¬ 
tion  approach. 

Channel  and  Algorithm  Sensitivities  —  There  are  two  important  conclusions  which  can 
be  drawn  from  the  numerical  penurbation  and  analytical  sensitivity  results.  First,  it  is  clear  from 
the  kernel  function  sensitivities  that  the  themial  wind  estimates  are  most  sensitive  to  uncertain¬ 
ties  in  the  kernel  functions  for  the  higher  wave  number  channels.  Therefore,  an  operational  im¬ 
plementation  of  the  optimal  estimation  algorithm  should  use  data  from  channels  sensing 
radiation  at  the  longer  wavelengths.  Of  course,  this  criterion  must  be  balanced  against  the 
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requiienient  to  use  a  sufFicient  number  of  channels  to  provide  sensitivity  to  the  horizontal  tem¬ 
perature  gradients  within  the  altitude  bands  oF  interest. 


Second,  analytical  results  showed  that  the  optimal  estimation  algorithm  was  much  less  sen¬ 
sitive  than  (3hring’s  method  to  kernel  function  uncertainties.  While  this  result  was  based  on  only 
one  case  (one  data  set  and  one  set  oF  channels),  it  is  consistent  w  ith  the  optimal  algorithm  design 
whicii  minimizes  a  Functional  (Appendi.x  Ci.  Eq.  3.1-19)  which  is  less  sensitive  to  the  kernel 
Function  than  that  used  by  Ohring  s  method  (Eq.  2.3-3). 


3.2  RECOMMENDATIONS 


The  following  specific  lecommendations  are  made  as  a  result  of  the  e.xperience  gained  in 
this  project.  The  ultimate  objective  toward  which  these  recommendations  lead  is  a  prototy'pe  of 
an  operational  implementation  of  the  optimal  estimation  algorithm  tor  themial  wind  profiles. 

i)  Complete  development  and  implementation  of  an  error  model  for  the  initial 
temperature  gradient  profile  uncertainty. 

ii)  Investigate  algorithm  sensitivity  to  realistic  temperature  profile  uncertainties. 

iii)  Use  supporting  satellite  images  to  evaluate  carefully  the  effect  of  clouds  in 
the  regions  surrounding  the  locations  at  which  the  thermal  wind  is  estimated 
in  a  representative  dttta  set. 

iv)  Conduct  a  statistical  evaluation  of  algorithm  perfomiance  over  a  suitable 
data  set. 
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APPENDIX  A 
RADIANCE  DATA 


This  appendix  contains  12  plots,  2  for  each  of  the  6  HIRS/2  channels  used  in  the  demon¬ 
stration.  Each  plot  shows  radiance  data  along  each  of  6  scan  lines,  either  the  upper  (more  north¬ 
erly)  or  lower  (more  southerly)  0  scans  spaced  approximately  40  km.  apart.  The  abscissa  (scan 
axis  location)  increases  in  the  easterly  direction.  Recall  that  the  scan  lines  are  rotated  approxi¬ 
mately  20“  to  the  north  of  due  east  (see  Fig.  2.1-2). 


RADIANCE  RADIANCE 


RADIANCE  RADIANCE 


Figure  A-3  Radiance  Data:  Channel  4-North 


Figure  A-4  Radiance  Data:  Channel  4-South 
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APPENDIX  B 

RADIANCE  GRADIENT  ESTIMATES 


The  six  plots  in  this  appendix  show  the  spatial  variability  of  the  radiance  gradient  vectors 
computed  from  the  data  using  a  planar  fit  over  a  7  x  7  grid  centered  at  each  of  the  25  locations 
in  each  plot.  Thus  an  adjacent  pair  of  gradient  vectors  (along  the  same  scan  line  or  the  same 
track)  were  computed  from  grids  containing  36  common  points.  The  plots  contain  qiialitciiive  in¬ 
formation  only;  magnitudes  were  normalized  so  that  the  longest  vector  fit  within  the  rectangle 
which  represents  the  (approximate)  40  km  x  40  km  grid.  The  unequal  scales  on  the  north  and 
east  axes  distort  the  true  angular  orientation  of  each  vector. 

Clearly,  spatial  variability  varies  significantly  between  channels.  Channel  3  gradients 
show  little  magnitude  variation  while  Channels  4.  5  and  6  have  uboiit  identical  large  magnitude 
variations.  Channel  15  varies  widely  in  both  magnitude  and  direction. 
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Figure  B-2  Radiance  Gradient  Vectors:  Channel  4 
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Figure  B-3  Radiance  Gradient  Vectors:  Channel  5 
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Fiyure  B-4  Radiance  (iradient  \  ectoi  s:  Channel  6 
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Figure  B-5  Radiance  Gradient  Vectors;  Channel  15 
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Figure  B-6  Radiance  Gradient  Vectors:  Channel  16 


APPENDIX  C 

KERNEL  FUNCTIONS 


The  two  plots  in  this  appendix  display  the  kernel  functions.  K,  (z),  used  in  the  demonstra¬ 
tion.  They  were  computed  from  the  formulas  presented  in  Appendix  G  using  the  Gibraltar 
measured  temperature  profile.  Fig.  2.2-1. 
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Figure  C-1  Kernel  Functions:  Channels  3.  4,  5.  6 
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APPENDIX  D 

RADIANCE  GRADIENTS  FROM  PLANAR  FITS 


The  Interim  Report,  Appendix  G.  presented  an  algorithm  for  computing  radiance  gradi¬ 
ents  from  the  array  of  measured  radiance  values  based  on  a  cubic  spline  fit  to  the  raw  data  field. 
Examining  results  of  that  algorithm  using  the  supplied  data  led  to  concern  regarding  the  vari¬ 
ability  of  the  results.  In  order  to  reduce  the  effects  of  measurement  errors  and  modeling  ap¬ 
proximations,  slopes  of  those  bicubic  spline  fits  were  averaged  over  specified  subsections  of  the 
complete  data  grid  (e.g.,  over  a  5  x  7  x  7  or  9  x  9  array).  Observed  variability  was  unexpect¬ 
edly  large  and  not  consistent  with  the  assumed  measurement  noise  models. 

To  avoid  further  loss  of  time  while  ttying  to  resolve  the  issues  raised  by  the  cubic-spline 
results,  an  alternative,  simpler  procedure  was  installed  based  on  planar  fits  to  the  observed  data 
values  over  specified  subgrids.  This  appendix  presents  the  equations  which  implemented  this 
procedure  and,  most  impottiintly.  the  corresponding  error  model  for  radiance  gradient  estima¬ 
tion  error  consistent  with  that  procedure.  Results  of  this  procedure,  summarized  in  Section  2. 
while  still  showing  some  \ariubility  with  sub-grid  size  and  location,  are  much  more  consistent 
with  the  assumed  radiance  measurement  error  model. 


D.l  PLANAR  FIT  ALGORITHM 

The  planar  fit  algorithm  is  based  on  the  following  problem  formulation; 

Given:  R(x,),  Radiance  measurements  at  one  wavelength  at  position 

X,.  i  =  1,...,N  ,  where  ^  =  (xj  .x:,)  define  a  two- 
dimensional  regular  grid 

Find  a  plane; 

f(;x)  =  a  +  bx,  +  cx2  =  d^  x'  (D.1-1) 
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where; 


(D.1-2) 


\'  =  (  1  x,  \2  )'  and  d'  =  (a  b  c) 

Such  that: 

=  X  |R(.Li  -  f'Csdl'  (D.1-3; 

i 

is  minimized. 


This  straightfonvard  least-squares  problem  has  tlie  well  known  solution 

d  =  .V‘  e 

in  which  the  3x3  matrix  A  and  the  vector  e  are  as  follows. 


e  =  ^  R(.x,).x.', 


(D.1-4) 


(D.1-5) 


(D.1-6) 


Equation  D.  1-1  defines  a  plane  whose  gradient,  (b  c).  is  an  estimate  of  the  radiance  gra¬ 
dient  at  the  center  of  the  data  grid  expressed  in  the  data  grid  coordinates.  These  coordinates 
were  along  the  track  and  scan  axes  (xi  and  x;  respectively)  of  the  satellite  data  collection  path. 
Since  the  wind  profile  estimation  will  be  conducted  in  an  (East.  Noah)  coordinate  system,  the 
above  gradient  vector  is  rotated  by  the  following  transformation: 


TCv--) 


(sin  »/'  cos  ip 
-  cos  ip  sin  ip 


(D.1-7) 


where  ip  =  counter-clockwise  angle  from  the  East  axis  to  the  scan  axis.  In  this  formulation,  it  is 
assumed  that  the  track  and  scan  axes  are  oahogonal.  The  actual  track  and  scan  axes  were 
slightly  non-onhogonal,  with  deviation  up  to  2°  in  some  parts  of  the  data  collection  region.  The 
orthogonal  approximation  will  cause  a  vei7  slight  eiTor  in  the  direction  and,  to  an  even  smaller 
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extent,  niagnittide  ot  calculated  radiance  gradients.  These  errors  are  far  below  the  noise  level  of 
the  radiance  gradient  estimates. 

Software  wa^  written  and  tested  implementing  Eqs.  D.1-4  -  D.1-7  for  user-selected 
square  subgrids  located  anywhere  within  the  complete  11x11  grid  selected  for  analysis.  Results 
are  presented  and  discussed  in  Section  2.1. 


D.2  PLANAR  FIT  ERROR  MODEL 


The  error  model  for  radiance  gradient  estimates  computed  by  the  above  procedure  is 
derived  from  two  sources,  the  measurement  error  model  for  the  radiance  data  and  the  modeling 
error  implied  by  the  planar  fit.  To  present  the  error  model  in  a  compact  fonn,  the  radiance  gra¬ 
dient  estimates  are  written 

VR'  =  HR  (D.:-l) 

here 

H  =  ;  T(v->  j  (X  X')-'  X  (D.2-2) 

R  =  R(ii)  R(:s:)  ...  R(:iN),'  (D.2-3) 

and 

=  2s’ 1  .s':  •  ■  S'n)  (D.3-4) 


Now, 


R  =  R,  +  R 


(D.2-5) 


where  the  subscript  T  indicates  the  vector  of  true  radiance  values  and  the  last  tenn  is  the  meas¬ 
urement  error  vector.  Therefore 

VR'  =  HRr  +  (D.2-6) 

and 

VR'  =  (VR'v,.  +  (D.2-7) 
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wliere  in  F.q.  D.2-7  the  gradient  is  expressed  as  a  siini  of  the  true  gradient  vector  plus  the  error 
term  whose  covtiriance  will  define  the  error  model.  If  the  true  radiance  values  actually  do,  in 
fact,  lie  in  a  plane,  then  there  is  no  modeling  error  (referred  to  earlier)  and  the  first  terms  on  the 
right  hand  sides  of  Eqs.  D.2-6  and  1  are  equal.  In  this  case,  the  error  model  is 

r  =  Cov  (VR)  =  H  C'ov  (R)H'  (D.2-8) 

Results  presented  in  this  repon  are  based  on  the  assumption  that  modeling  error  can  be  ne¬ 
glected  (so  that  Eq.  D.2-S  holds). 

The  covariance  of  the  radiance  data  vector  error.  R.  in  this  eLjuation  is  exactly  the  same 
as  the  matrix  discussed  in  Section  2  of  the  Interim  repon.  All  analysis  results  reponed  in  this 
document  used  a  <t>  matrix  based  on  parameter  values  (measurement  error  variances  and  corre¬ 
lation  distances)  stated  in  that  report. 
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APPENDIX  E 

PROOF  OF  ALGORITHM  VALIDITY 


The  objective  of  this  proof  is  to  demonstrate  the  consistency  of  wind  profile  estimates, 
i.e.,  that  in  the  absence  of  noise  and  numerical  error,  the  algorithm  will  produce  the  correct 
thermal  wind  estimates.  A  special  case  of  consistency  was  shown  earlier  in  the  Interim  Repon 
(Appendix  G).  Though  the  simple  lineai'  form  assumed  for  the  temperature  field  in  that  special 
case  was  useful  in  the  simulation-based  check  of  the  wind  profile  estimation  algorithm,  it  was. 
in  fact,  not  required  in  the  analytical  proof  that  wind  estimates  converge  to  the  true  thermal  wind 
in  the  limit  as  the  radiance  gradient  error  covariance  matri.v  T  approaches  the  zero  matrix. 

The  critical  step  in  the  proof  is  recognition  that,  for  any  continuously  differentiable  tem¬ 
perature  field.  T(y.z)  (since  we  are  presently  examining  only  the  east  component  of  the  thermal 
wind,  the  x  dependence  is  not  indicated),  the  following  relationship  holds; 

g  =  ^  (E-1) 

where  both  sides  of  Eq.  E-1  were  defined  in  Section  2.2.1.  Note  that  distinct  symbols  are  re¬ 
quired  on  the  two  sides  of  Eq.  E-1  because  g  is  based  on  radiance  measurements  whereas  9  is 
based  on  the  assumed  temperature  gradient. 


To  prove  that 

lim  u  =  u” 

r  ^  0 

it  is  sufficient,  via  Eq.  2.2-7  and  2.2-8,  to  prove  that 

lim  g*  =  1 

r  ->1) 

I.e..  that 

lim  o'  P-'q  =  1 
\  -*  0 

by  Eq.  E-1 . 


(E-2) 


(E-3) 


(E-4) 
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The  remainder  of  the  proof  follows  verbatim  that  for  the  linearly  varying  temperature 
field.  Applying  the  Shemian-Morrison  formula  (Ref.  3)  to  Eq.  2.2-4.  P  '  is  given  by 


p-i  ^  p-i 


( _ ' _ 


r* 


(E-5; 


Now.  the  left  side  of  Eq.  E-4  can  be  e.xpressed  as 


lim  {0^  P-‘  ^1)  = 

r  -*() 


lim 

r  0 


J 


r  f)'  r 


which  becomes 


lim 

r  -  0 


(O'  P  ' 


=  1 


(E-6) 


fE-7) 


Thus  the  thermal  wind  estimate  u  approaches  the  true  value  u"  as  the  radiance  gradient 
eiTor  covariance  matri.x  approaches  zero.  This  result  not  only  provides  theoretical  justification 
for  the  algorithm,  but  also  enables  a  check  for  vertical  discretization  and  numerical  integration 
errors. 
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APPENDIX  F 

SENSITIVITY  EQUATIONS 


Although  the  most  important  results  regarding  robustness  of  the  new  optimal  estimation 
algorithm  are  based  on  actual  algorithm  executions  using  perturbed  parameter  values,  sensitiv¬ 
ity  equations  based  on  linearized  equations  provide  imponant  theoretical  results  in  a  cost- 
efficient  manner.  These  linearized  sensitivity  equations  are  presented  in  this  appendix  for  both 
the  optimal  estimation  approach  and  for  Ohring’s  method.  Numerical  results  based  on  these 
equations  are  presented  in  Section  2.4  aloi.g  with  algorithm  penurbation  results. 


F.l  SENSITIVITY  OF  THE  OPTIMAL  ESTIMATION  ALGORITHM 

Sensitivity  equations  can  be  derived  for  any  of  the  algorithm  parameters  or  data  values. 
This  section  presents  equations  only  for  the  sensitivity  with  respect  to  values  of  the  kernel  func¬ 
tions,  Ki(z),  because  those  functions  depend  on  the  most  imponant  input  information,  namely 
the  nominal  temperature  profile,  the  Planck  functions,  and  the  absorption  profile  (see  Appen¬ 
dix  G).  Sensitivity  equations  with  respect  to  the  assumed  temperature  gradient  profile  could  be 
derived  in  a  similar  manner. 

There  is  one  imponant  technical  point  in  the  following  derivation  that  should  be  noted. 
The  kernels,  Kj(z)  are  functions  of  z,  so  a  functional  analysis  viewpoint  (using  a  Frechet 
derivative,  for  example)  should  be  used  to  express  sensitivity  with  respect  to  the  “kernel  func- 
tions".  As  an  alternative,  simpler  approach,  the  following  derivation  is  based  on  an  approxi¬ 
mately  uniform  fractional  penurbation  throughout  jO.  zj). 

The  sensitivity  equations  are  based  on  a  logarithmic  derivative  of  the  estimates,  u. 
(Eq.  2.2-2)  with  respect  to  each  of  the  kernel  function  components.  So.  there  is  a  n  x  m  matrix 
of  sensitivity  functions  defined  by: 


The  element  of  S  is  given  by 


S„  =  - 


R  Kj  ()  c(Z|)’ 
f  ui  (?i)  <'>k, 


i  =  1, 
j  =  1- 


....  n 
....  m 


(F.1-2) 


By  e.xtensive  bur  straightfonvard  manipulation,  the  following  formula  for  the  complete 
sensitivity  matri.x  results: 

R 


S  =  a  I  Uj)  ^t2  t|  -p  (  j  K|> 


where 


a,  =  [  '  ---  {C)dC 

Jo  fly 

U()  =  Diag  i(l/u(7,)l 

t|  =  P'’  g  .  ti  =  (i  -  f  fj 
(C  is  given  by  Eq.  2.2-3  and  a  by  Eq.  2.2-6), 

/J  =  t, 


and 


Ki)  =  Diag(kj) 


where  the  diagonal  elements  are  constants. 

kj  =  —  f  '  I  Kj  (0  I  dC 

Z I  J  n 


(F.1-3) 

(F.1-4) 

(F.1-5) 

(F.1-6) 

(F.1-7) 

(F.1-8) 

(F.1-9) 


F.2  SENSITIVITY  OF  OHRING’S  METHOD 

The  approach  used  to  develop  a  sensitivity  equation  for  Ohring’s  method  parallels  that  in 
the  previous  section,  with  a  few  necessary  modifications.  Since  Ohring’s  method  is  formulated 
to  yield  layer  estimates  (Eq.  2.3-1)  we  define  the  vector  of  summed  estimates  analogous  to  that 
produced  by  the  optimal  estimation  approach  as  follows: 

MO  =  L  uj.  (F.2-1) 


F-2 


where  U],  is  the  vector  of  estimates  calculated  using  Eq.  2.3-1,  and  the  matrix  L  has  elements 


ij 


( 


^  J 

<  j 


(F.2-2) 


So,  the  sensitivity  matrix  for  Ohring’s  method  is 

S  =  -(R/0  U(,  L  (tiotjo  -  A)  Co)  Ko  (F.2-3) 


where  Uois  a  diagonal  matrix  analogous  to  Eq.  F.1-5  but  containing  Ohring’s  summed  esti¬ 
mates,  Uio,  and  Kd  is  again  given  by  Eqs.  F.1-8  and  F.1-9.  The  other  terms  are: 

tio  =  ^  -  Co  k  (F.2-4) 


where 


Az,  =  z,  +.  1  -  z,.  i  =  1 . n. 


k 


K(C)t]^, 


t20  —  M"'  Ry 


where  M  is  given  by  Eq.  2.3-6, 


and  finally, 


(F.2-5) 


(F.2-6) 


(F.2-7) 


(F.2-8) 


(F.2-9) 


where  Cj ,  i  =  1 .  n,  are  given  by  Eq.  2.3-5. 

The  similar  forms  of  the  Ohring  and  optimal  estimate  sensitivity  matrices  should  not  be 
surprising  because  of  the  outwardly  similar  form  of  the  estimators.  That  their  numerical  values 
are  much  different  is  shown  in  Section  2.4. 
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Research  efforts  under  the  referenced  con¬ 
tract  have  focused  on  developing  a  practical  ap¬ 
proach  for  estimating  thermal  wind  profiles  di¬ 
rectly  from  satellite  radiance  data.  The  approach, 
if  successful,  can  be  adapted  for  operational  use 
with  the  addition  of  at  least  one  single  level  wind 
measurement  to  yield  an  estimate  of  the  vertical 
wind  profile.  In  this  interim  report  the  baseline 
algorithm,  incorporating  a  statistical  model  of 
satellite  radiance  gradient  errors,  is  reviewed  and 
preliminary  simulation  results  are  presented.  A 
final  report,  due  in  December  1988,  will  provide 
the  final  algorithm  and  demonstration  results. 
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INTRODUCTION 


1  . 


The  need  for  determining  wind  fields  in  data-sparse 
regions  has  prompted  several  investigations  into  the  use  of 
satellite  radiance  data  for  wind  profile  determination.  Ref¬ 
erences  1,  2,  and  3  describe  methods  which  compute  thermal 
wind  directly  from  radiance  measurements,  avoiding  the  error- 
prone  path  of  first  determining  temperature  profiles  (the 
inversion  problem)  and  then  computing  thermal  wind  via  the 
thermal  wind  equation.  The  present  effort  extends  previous 
work,  (notably  Ref.  1)  by  formulating  an  optimal  estimation 
algorithm  based  on  a  statistical  error  minimization  criterion. 
Statistical  models  of  error  sources  in  the  data  collection  and 
data  analysis  processes  are  an  integral  part  of  the  formula¬ 
tion.  In  this  report,  only  the  radiance  gradient  error  model 
has  been  defined.  The  algorithm  structure  can  accommodate 
error  models  for  temperature  gradients  and  f requency-dependent 
weighting  functions  as  well.  Descriptions  of  these  models 
will  appear  in  the  final  report. 


Background,  motivation,  and  discussions  of  previous 
work  related  to  this  effort  are  found  in  Ref.  4.  The  work 
tasks,  as  outlined  in  Ref.  4,  are  listed  in  Figure  1-1.  Pro¬ 
gress  achieved  to-date  for  each  task  is  indicated  by  shaded 
bars  . 


The  remainder  of  this  report  is  organized  as  follows. 
Section  2  provides  a  description  of  the  satellite  weighting 
function  models  used  in  this  effort.  A  description  of  the  ra¬ 
diance  gradient  error  model  is  also  provided.  Section  3  de¬ 
tails  the  baseline  algorithm  formulation  including  all 
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Figure  1-1  Progress  in  Development  of  the  Statistical 
Wind  Profile  Estimator 

assumptions  and  equations.  An  extension  to  the  baseline  algo¬ 
rithm,  based  on  a  closed-form  iterative  approach,  is  described 
as  well.  In  Section  4  the  data  to  be  used  in  producing  pre¬ 
liminary  demonstration  results  are  identified  aiid  described, 
and  simulation  results  for  algorithm  and  software  verification 
are  presented  and  discussed  in  Section  5.  A  summary  is  pro¬ 
vided  in  Section  6. 


For  use  during  algorithm  development  and  testing, 
parametric  forms  for  the  baseline  kernel  functions  and  the  ra¬ 
diance  gradient  error  model  have  been  defined.  In  addition  to 
their  utility  during  algorithm  checkout  and  demonstration, 
these  models  provide  the  means  to  determine  estimate  sensi¬ 
tivities  to  various  model  components.  The  parametric  models 
for  the  kernel  functions  and  the  radiance  gradient  error  model 
are  described  in  .Sections  2.L  and  2.2,  respectively. 

2.1  SPECIFICATION  OF  THE  KERNEL  FUNCTIONS 

AS  will  be  discussed  in  Section  3  the  theoretical  ra¬ 
diance  measured  by  the  satellite  radiometer  in  each  specific 
channel  is  computed  from  the  radiative  transfer  equation  in¬ 
volving  a  vertical  integral  of  a  product  of  temperature  and  a 
f requency-dependent  kernel  function.  The  kernel  or  weighting 

function  K.  (z)  associated  with  the  i^*^  channel  of  the  satellite 
1 

radiometer  is  given  in  Ref.  4  as 

dB;  dT; 

Ki(z)  ^  [T^  (2.1-1) 

where  B.  is  the  Planck  function  at  the  i^*^  frequency  or  wave 
number,  T  is  temperature,  z  is  an  altitude  variable,  T^  (z)  is 
a  given  nominal  temperature  profile,  and  t.  is  the  transmit¬ 
tance  function  at  the  i^^  frequency. 


The  Planck  function  for  a  discrete  wave  number  v.  is 

1 

given  as 


'-1 

B.  (T)  =  -  (2.1-2) 

exp  (c^  V.  /T)  -  1 

where  B  is  black  body  radiance  in  W  m-2  sr'’-  cm,  T  is  absolute 
temperature  (deg  K)  .  and  c^  and  c^  are  constants  (c^^  =1 .  191 1  10-8 
W  m-2  sr-i  cm^  and  C2=l.439  cm  deg).  Differentiating  Eq.  2.1-2 
with  respect  to  temperature  yields 

HB  exp  (c  V.  /T) 

_2£i  =  _ i _ ± _ I _ LJ _  (2.1-3) 

T2[exp(C2V.  /T)-l 


Values  for  this  term  at  each  wave  number  of  interest  are  com¬ 
puted  using  the  nominal  temperature  profile. 

The  transmittance  t, (z)  at  wave  number  v.  from  alti- 

1  1 

tude  z  to  the  "top"  of  the  atmosphere  is  a  complex  function  of 
atmospheric  constituents.  In  this  effort  the  vertical  deriva¬ 
tives  of  T.  needed  to  evaluate  the  kernel  functions  (Eq.  2.1-1) 
are  represented  by  analytic  expressions  fit  to  published 
curves.  Plots  of  the  derivative  function 


dfi 


(2.1-4) 


for  the  HIRS/2  are  contained  in  Figures  2.1-1  and  2.1-2  (see 
Ref.  5).  Parameter  values  a,  b,  c  at  each  v,  for  which  the 
following  functional  form 

4)(p)  =  a  p^^  exp(-cp)  (2.1-5) 
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Figure  2,1-1  Vertical  Derivatives  of  Transmittance 

Functions,  Channels  1-8  ^ ^ 


dr  'd(p"') 

Figure  ?  1-2  Vertical  Derivatives  of  Transmittance 

Funetirins,  Channels  13-19 


j-  •) 


best  fits  Eq.  2.1-4  have  been  determined.  The  best  fit  crite¬ 
rion  used  is  the  quadratic  function 


n 

^[t:(p)-4)(p,)]'^  (2.1-6J 

1  j  j 

j  =  l 

where  the  n  pressure  levels  are  preselected.  Sample  fits  are 
illustrated  in  Figure  2.1-3.  As  a  rule,  the  fit  is  best  for 
roughly  bell-shaped  curves,  e.g.  those  corresponding  to  chan¬ 
nels  2  and  15.  The  fit  for  the  curve  corresponding  to  channel 
6  is  improved  by  weighting  points  in  the  250-1000  mb  range 
more  heavily,  as  illustrated  by  the  dotted  line.  From  these 
parametric  forms,  it  follows  that 

=  L  ,  f"'-  =  f-f  exp(-cp)  (2.1-7) 

dz  d ( In  p)  yi  J 

The  product  of  Eqs .  2.1-3  and  2.1-7  is,  for  each  frequency, 
the  desired  kernel  function. 


2.2  RADIANCE  GRADIENT  ERROR  MODEL 

In  this  section  a  matrix  expression  for  the  fre¬ 
quency-dependent  radiance  gradient  error  model  is  derived. 
Publisheu  values  of  HIRS/2  radiance  measurement  error  standard 
deviations  are  directly  incorporated  using  a  prescribed  ana- 
l..tic  measurement  error  covariance  model.  The  geometry  of  the 
gradient  computation  is  contained  in  the  gradient  error  model 
by  modeling  the  specific  bicubic  spl’ne  derivative  operations 
made  on  the  radiance  data.  For  each  frequency,  radiance 
measurements  form  a  sea  1 ar - va 1 ucd  function  r  defined  on  a  pla¬ 
nar,  N  by  N  grid  as  pictured  in  Figure  2.2-1.  A  standard 

X  y  ° 
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Figure  2.1-3  Vertical  uerivatives  of  Transmittance 

Functions  (Solid  Lines)  and  Parametric 
Fits  (Broken  Lines) 
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Figure  2.2-1  Radiance  Grid 

bicubic  spline  routine  yields  radiance  gradients,  for  each 
frequency,  at  any  preselected  grid  point.  It  is  shown  in  Ap¬ 
pendix  A  that  the  radiance  gradient  at  grid  point  (j,k.)  is  a 
linear  transformation  of  the  vector  jr  of  all  radiance  measure 
ments  over  the  grid,  i.e. 


=  F-  1,  r 

^  j .  k  — 


(2.2-1) 


where 


^  =  <''1,1  '■1.2 


r  r  r  r 

1 .  .'r,  2.1  2.2  2 .  N  ^ 


r  r 

•Ny,  1  .\y.2  •  ■  ■  Sy.N.^ 

(2.2-2) 


and  F.  ,  is  a  certain  2  by  (N  *  N  )  matrix  which,  as  indicated 

j  ^  Ic  X  y 

notat ional ly ,  is  dependent  on  the  gridpoint.  Equation  2.2-1 
provides  an  analytically  convenient  form  for  deriving  the 


radiance  gradient  error  model  given  a  radiance  measurement  er¬ 
ror  model,  especially  in  the  case  in  which  spatial  correla¬ 
tions  are  modeled.  Since  Eq .  2.2-1  yields  the  gradient  at  one 
point,  an  averaging  procedure  has  been  defined  to  reduce  data 
errors  and  smooth  results,  if  necessary.  The  average  radiance 
gradient  over  a  prespecified  subgrid  of  points  is  given  by 


(2.2-3) 


=  F  L 


where  the  summation  extends  over  the  subgrid  and  where 


j  k 


(2.2-4) 


The  averaging  process  yields  a  more  representative  radiance 
gradient  in  the  subgrid  region.  It  follows  that  the  radiance 
gradient  error  covariance  matrix  F,  defined  by 


F  =  E  -I  £r  vr 

I  -“-^avg 


(2 . 2-5  ) 


where  ^ave  is  the  error  in  the  computed  gradient,  can  be 
written  as 


F  =  F  4)  F 


(2 . 2-6) 


where  O  denotes  the  radiance  measurement  error  covariance  ma¬ 
trix.  The  remainder  of  this  section  is  hence  devoted  to  the 
specification  of  <P . 


Under  the  reasonable  assumptions  that  radiance  meas¬ 
urement  errors  for  one  frequency  are  isotropic  in  space  and 
possess  an  autocorre 1  at i on  function  that  decays  exponentially 
with  horizontal  distance,  it  follows  that 


where 

r  =  r  -  r  (2 . 2-8) 

the  difference  between  measurement  and  truth,  a  and  b  refer  to 
spatial  locations,  d  denotes  the  Euclidean  distance  between 
a  and  b,  D  is  a  correlation  distance  parameter  and  o  is  the 
instrument  noise  variance  for  the  channel.  For  preliminary 
algorithm  work,  the  parameter  D  is  taken  to  be  300  km.  Error 
standard  deviations  from  Ref.  6  are  listed  in  Table  2.2-1. 
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TABLE  2.2-1 

HIRS/2  RADIANCE  ERROR  STANDARD  DEVIATIONS 


CHANNEL 

CHANNEL  WAVE  NUMBER 

ERROR  STD.  DEV. 

(cm“^  ) 

(mW/m^  sr  cm"^  ) 

1 

669 

3 . 00 

2 

680 

0 . 67 

3 

690 

0 . 50 

4 

7  03 

0.31 

5 

716 

0.21 

6 

733 

0 . 24 

7 

749 

0 . 20 

8 

900 

0  .  10 

13 

2190 

0 . 006 

14 

2210 

0 . 003 

15 

2240 

0 . 004 

16 

2270 

0 . 002 

17 

2420 

0 . 002 

18 

2515 

0 . 002 

19 

2660 

0.001 
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3.1  THE  BASELINE  ALGORITHM 

The  radiative  transfer  equation  and  the  thermal  wind 
equation  are  the  key  elements  in  the  derivation  of  the 
baseline  algorithm.  The  radiative  transfer  equation  relates 
the  radiance,  R.  ,  measured  by  a  sate  1 1 i te - borne  radiometer  at 
wave  number  v.  ,  to  the  atmospheric  temperature  profile.  The 
following  simple  form  of  the  radiative  transfer  equation, 
based  on  a  given  temperature  profile  (e.g.,  from  a  forecast 
field),  is  derived  in  Ref.  4: 

Ri  =  Rf  +  ^  Ki  (z)  T(x,  y,  z)  dz  (3.1-1) 

where , 

x,y  define  zonal  and  meridional  directions,  respectively 

R°  is  a  constant 
1 

z.^  defines  the  'top  of  the  atmosphere' 

K. (z)  is  an  approximately  known  weighting  function  at  wave 

number  v.  (see  Section  2.1) 

T(x.y,z)  is  the  temperature  profile  at  location  (x,y) 

z  is  the  vertical  coordinate  defined  by 

z  =  ln(po/p) 
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where  p  denotes  pressure  and  p^  denotes  surface  pressure. 

Assuming  that  horizontal  variations  of  are  negli¬ 
gible  (following  the  treatment  given  in  Ref.  4),  the  radiance 
gradient  is  determined  from  Eq .  3.1-1  as 

VRi  (x,y)  =  ^  K,  (z)  VT  (x,y,  z)  dz  (3.1-2) 

Equation  3.1-2  represents  a  two-vector  with  x  and  y  components 


^  =  r  K, 

(z) 

aT 

(z)dz 

(3. 1-3) 

d.x  Jo 

ax 

^  =  r  K, 

(z) 

aT 

(z)dz 

(3.1-4) 

dy  Jo 

ay 

where  the  dependence  on  horizontal  position  has  been  sup¬ 
pressed  . 


The  thermal  wind  equation  relates  the  thermal  wind, 

X-j- .  i  .e.,  the  vector  difference  between  geostrophic  winds  (Y^^) 
at  two  levels,  to  the  horizontal  temperature  gradient 

V-r  (Z|  zj  =  Vf;  (zj  -  V,i  (i^)  =  y  f  “(k  X  VT)  dz  (3.1-5) 

t  JZf 

whe  re 

R  is  the  gas  constant 

f  =  2f2Sin(4))  is  the  Coriolis  constant  at  latitude  4)  (Q 
is  the  earth  rotation  rate) 
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k  is  a  unit  vector  in  the  vertical  direction. 

The  horizontal  gradient  in  this  expression  is  within 
constant  pressure  surfaces  as  opposed  to  a  fixed  geometric  al¬ 
titude.  As  discussed  in  Ref.  4,  however,  the  difference  is 
assumed  negligible  in  this  effort.  The  scalar  components  of 
V  are  given  by 


Ut(z^,zJ  =  Ur,(Zu)  -  Uo(z^) 


f  dy 


(3.1-6) 


VT(Zi.Zu)  =  Vo(Zuj  -  V(;(Zj,) 


R 

f 


It  d\ 


dz 


(3.1-7) 


where  u  and  v  denote  the  zonal  and  meridional  components,  re¬ 
spective  ly  . 

The  possibility  of  estimating  thermal  wind  d i rect Iv 
from  radiance  data  is  suggested  by  the  similarity  between 
Eqs .  3.1-2  and  3.1-5.  Both  involve  vertical  integrals  of  the 
temperature  gradient  profile.  However,  the  integral  in 
Eq .  3.1-2  is  over  the  entire  atmosphere  and  is  weighted  by  K.  , 
while  that  in  Eq .  3.1-5  is  unweighted  and  limited  to  the  layer 
defined  by  ( zj  ,  )  .  In  light  of  this,  a  function  is  defined 


W(z) 


1  Zjj  <  Z  <  Zu 

0  otherwise 


(3.1-8) 


For  rotational  convenience  the  dependence  on  the  interval 
( zj  . )  is  suppressed  throughout  the  remainder  of  the  report. 
Equation  3.1-5  can  now  be  written  as 
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-  T  W  (z)  [  k  X  ^  ]  dz 
f  Jo 


(3.1-9) 


The  zonal  and  meridional  scalar  components  are  given  respec¬ 
tively  by 


Uj 


R 

f 


■ir 

.  o 


W(z)  dz 
by 


(3.1-10) 


Vt  = 


(3.1-11) 


Now,  a  linear  estimator  is  defined  using  computed  (that  is, 
estimated)  radiance  gradients  for  available  frequency  channels 


Ut 


R 

T 


I  . 

•  -7-  C, 

I  by 


(3.  1-12) 


where  the  c^^  are  constants  to  be  optimally  determined.  Hence¬ 
forth,  for  simplicity,  only  the  zonal  component  will  be 
treated  here.  An  analogous  procedure  for  the  meridional  com¬ 
ponent  is  evident  using  radiance  derivatives  in  the  x  direc¬ 
tion.  The  computed  radiance  gradients  consist  of  true  plus 
error  quantities  (~  denotes  an  error  quantity) 


dRj  r(Ri  (9R| 

by  by  by 


(3.1-13) 


For  convenience  in  the  algebraic  manipulations  which 
follow,  the  radiance  gradient  terms  and  kernel  functions  for 
all  available  frequencies  are  organized  into  vectors,  i.e.. 
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g  = 

dy 

dy 

:  - 

K  = 


Thus  Eq.  3.1-12  is  expressed  in  the  form 


»  R  -  j 

Ut  =  -J  g  c 


where 


(3.1-14) 


(3.1-15) 


(3.1-16) 


The  statistical  estimation  problem  consists  of  determining  the 
'' oef f icients  q,  to  minimize  the  estimation  error.  For  the  zo¬ 
nal  component  this  is  gi .  ..  by  the  difference  between 
Eqs .  3.1-15  and  3.1-10 


(3.1-17) 


Using  Eqs.  3.1-4  and  3.1-13,  Eq .  3.1-17  is  expressed  as 

Ut  =  -  Y  c  -  f  ^  [W  (z)  -  ^  dz  1 


(3. 1-18) 


The  statistical  minimization  criterion  is  represented  by 


J  =  E  {u\} 


(3 . 1-19) 


and  is  quadratic  in  By  differentiating  Eq.  3,1-19  with  re¬ 

spect  to  Q.  and  setting  the  result  equal  to  zero,  the  minimiz¬ 
ing  solution  is  found  to  be 


c  =  P  ’  q 


(3  .1-20) 


where 

p  =  r  +  Of 


(3.1-21) 


q  =  aO 


(3 . 1-22) 


In  Eq .  3,1-21,  the  F  matrix  represents  the  covariance  of  the 
1  \diance  gradient  errors  as  described  in  Section  2.2.  The 
vector  £  and  the  scalar  a  involve  the  temperature  gradient  and 
the  satellite  weighting  functions  and  are  computed  by  the  fi..l- 
lowing  numerical  quadratures; 


(3.1-23) 


(3  .1-24) 
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where  is  the  number  of  vertical  levels  and  the  symbol  "=" 
denotes  approximate  equality.  The  root  mean- square  -  error 
(rmse)  in  the  thermal  wind  estimate  is  given  by 


rmse 


(  2 

f 


P- 


3)’ 


(3.1-25) 


It  is  worth  noting  that  the  minimization  criterion 
results  in  a  biased  thermal  wind  estimate.  This  is  seen  by 
taking  the  expected  value  of  Eq.  3.1-18.  The  bias  is  found  to 
be 

(utJb  =  V  f  ^  t  W(z)  -  K‘'c]~dz  (3.1-26) 

f  Jo  3y 


.\lthough  an  unbiased  estimate  could  be  determined,  it  would 
not  be  optimal  in  the  sense  of  minimum  mean  -  square  -  error . 


3.2  ITERATIONS  AND  STEADY-STATE 

Reference  4  suggested  a  method  for  extending  the  al¬ 
gorithm  in  Section  3.1.  As  a  consequence  of  Eq .  3.1-6,  it 
makes  sense  to  update  the  guess  of  the  temperature  derivative 
profile  by 


£ 

R 


dz^, 


u-r  (z) 


g'  P-'  B 


(3.2-1) 

(3.2-2) 


where  the  second  equality  follows  from  Eqs .  3.1-15  and  3.1-20. 
The  new  guess  of  the  temperature  derivative  profile  could 
then,  in  turn,  be  used  to  compute  a  new  thermal  wind  by 


O-iy 


Eq.  3.1-12  and  hence  another  update  of  the  temperature  deriva¬ 
tive  by  Eq.  3.2-2.  This  iterative  process  could  be  carried 
out  indefinitely,  assuming  it  converges  eventually  to  a  physi¬ 
cally  meaningful  function.  The  validity  of  this  assumption 
has  been  the  focal  point  of  a  thorough  investigation  of  the 
iterative  approach,  of  which  only  conclusions  are  reported 
here  . 


The  setting  defined  in  Section  3.1  (considering  only 
the  zonal  wind  component)  is  a  special  case  of  the  setting  de 
fined  in  Ref.  4  (considering  both  zonal  and  meridonal  wind 
components  s imviltaneous Iv  and  minimizing  a  sum  of  squares  cri 
terion) .  For  the  sake  of  maximum  generality,  we  discuss  it¬ 
erative  investigation  conclusions  in  the  latter  setting.  To 
define  circumstances  under  which  the  iterative  procedure 
yields  meaningful  results  requires  changing  the  notation  in 
Section  3.1  to  fit  the  more  general  vector  case.  Equa¬ 
tion  3.1-15  becomes 


*  R 

Vj  =  --  G  c 

_T  j 


G  = 


:  a  2 

X  n  matrix 

0  f 

dR] 

CIR2 

dy 

by 

c)Ri 

r)R2 

r3x 

i 

i  ven 

by  Eq .  3  . 

1  -  2 

r  + 

0,  0,' 

P 

q  =  0  a 


(3.2-3) 


(3.2-4) 


(3.2-5) 

(3.2-6) 
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and  0  and  a  are  the  matrix  and  vector  generalizations,  respec¬ 
tively  of  a  defined  by  Eqs .  3.1-23  and  3.1-24.  The  itera¬ 
tion  given  by  Eq .  3,2-2  is  now  replaced  by 

VLm  (z)  =  G  P-'  0i  (z)  (3.2-7) 

where  i  is  the  iteration  index. 

Many  approaches  exist  in  the  mathematical  literature 
for  determining  conditions  necessary  and/or  sufficient  for  a 
function  sequence 

(ZLiz)  }  (3.2-8) 

to  converge.  In  this  case,  the  most  promising  approach  in¬ 
volves  looking  not  directly  at  the  above  function  sequence, 
but  instead  at  an  associated  nxn  matrix  sequence 

(3.2-9) 

where 

Xi  =  0i  0i‘^  (3.2-10) 

Note  that  each  X,  has  rank  at  most  2  since  0.  is  nx2 .  The 

1  1 

advantages  of  a  finite-dimensional  linear  algebraic  formula¬ 
tion  over  an  infinite-dimensional  functional  analytic  formula¬ 
tion  will  become  obvious  as  the  discussion  progresses.  Since 
combining  Eqs.  3.1-23  and  3.2-7  yields 

0i‘;,  =  G(r  +  0,0^  ' .  (3.2-11) 

it  follows  that 
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(3.2-12) 


X,.,  =  Xi  (r  +  X,)-'  G^G  (r  +  Xi)-'  x,  , 


i.e.,  the  matrix  recursion  involves  no  integrations.  Note 
that  if 


1VI,>  VI 


(3.2-13) 


then 

X  =  lim  Xi 
i  -►  » 


(3 .2-14) 


uniform  con- 
VT  must  sat- 

n:  =  G  ( r  +  X)-’  x(r  +  x)-'  G"^  VY  (3.2-15) 


exists,  i.e.  ,  convergence  of  {X,}  is  necessary  for 
vergence  of  {VTj};  moreover,  the  limiting  function 
is  fy 


i.e. ,  ^  is  a  fixed  point  under  the  action  of  a  certain  2x2 
matrix  based  on  X.  In  particular,  the  direction  of  VT  is  con- 


(with  respect  to 


unless 


G  (r  +  x)  ’  x(r  +  x)-‘  G^  =  I  (3.2-16) 

No  analogous  restrictions  on  the  length  of  rr  have  yet  been 
found  (approximate  schemes  are  under  consideration).  The  im¬ 
portance  of  obtaining  the  limit  X  (in  closed  form,  if  possi¬ 
ble)  is  nevertheless  evident. 

The  following  Theorem  provides  conditions  sufficient 
for  (Xj)  to  converge. 
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Theorem  Suppose  that  there  exist  n-vectors 
a.  b,  u,  y 


and  scalars 


Xo  >  0,  Vo  >  0 

i 

satisfying  the  two  decomoosab i 1 i tv  conditions 


G^G  =  a  +  b 

Xo  =  Xo  uu^  +  Vo  yy‘ 

and  the  two  orthogonality  conditions 

y^  r*’  u  =0 
0(3  +  /Sy  =  0 

where 


(3.2-17) 


(3. 2-18) 


a  =  r*'  a 

y  =  u'f  r’  b 


=  y*'  r‘  b 
(5  =  r’  a 


(3.2-18) 


Then 


X  =  lim  Xj 

i  -►  ■» 

exists.  Moreover,  there  exist  x  >  0  .>..d  y  >.  0  such  that 
X  =  xuu^  +  yyy' 


Explicit  formulas  can  be  written  for  x  and  y  under  each  of 
four  cases . 

For  convenience,  let 
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A  =  r’  u 

=  v’'  r'  y 

Formulas  are  written  out  for  x  in  each  case.  Similar 
for  y  are  obtained  by  replacing  a,  y.  A,  Xq  by  <5,  ju,  yo 
tively  in  each  expression 

Case  i )  If 

0  <  4A  <  +  '/ 

and 

+  '/  -  2A  -  yja^  +  -  4A  <  2A^Xo, 

then 

-  2A  +  v'rt^  +  y^  +  y^  -  4A 

X  = 


Case  i  i ;  If 

0  <  4A  <  +  ’/ 

and 

a’  +  '/  -  2A  -  v'a^  +  T  +  y^  -  4A  =  2A^Xo, 

then 

X  =  Xo  . 


Case  i  i  i )  If 

0  <  4A  =  +  y^ 


and 

+  y^  -  2A  <  2A’xo  , 

then 


X 


+  y^  -  2A 


Case  iv)  Otherwise  x  =  0. 


formulas 

respec- 
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Thus  one  observes  from  the  Theorem  an  "upper  bound" 
on  the  covariance  matrix  P  (in  terms  of  the  two  quantities  A 
and  n)  for  which  nonzero  limits  x  exist.  This  upper  bound  may 
be  thought  of  as  a  "reasonableness"  bound,  beyond  which  the 
iteration  procedure  converges  to  physically  unrealistic  pro¬ 
files. 

We  wish  to  emphasize  the  importance  of  closed  form. 
solutions  as  presented  in  the  Theorem  statement.  The  itera¬ 
tive  procedure  for  estimating  temperature  gradient  profiles 
need  not  be  iterative  if  such  formulas  are  always  applicable. 
The  advantage  of  such  formulas  is  that  they  would  not  be  sub¬ 
ject  to  possible  slow  convergence  rates  which  have  frequently 
occurred  during  computer  experimentation  using  artificial  in¬ 
puts  (and  which  might  be  hazardous  and  expensive  if  programmed 
into  a  wind  profile  estimation  procedure). 


4.  SUPPORTING  DATA  FOR  ALGORITHM  DEMONSTRATION 


The  radiance  data,  rawinsonde  data,  and  temperature 
field  data  described  in  this  section  have  been  acquired  for 
purposes  of  algorithm  evaluation  and  demonstration.  No  data 
processing  results  have  been  prepared  for  this  interim  report, 
but  demonstration  results  will  be  featured  in  the  final  report 
due  at  the  end  of  the  contract  period. 


4.1  RADIANCE  DATA 

Radiance  data  from  the  High  Resolution  Infrared  Ra¬ 
diation  Sounder/2  (HIRS/2)  of  the  NOAA-7  polar  orbiting  satel¬ 
lite  has  been  obtained  for  algorithm  demonstration  purposes. 
The  device  is  a  part  of  the  TIROS  Operational  Vertical  Sounder 
(TOVS)  instrument  package  and  senses  radiation  in  20  channels 
in  the  infrared  spectrum.  The  data  set  covers  the  region  10®W 
-  5°E  and  SO^N  -  43®N  (Western  Europe  and  a  small  part  of 
North  Africa)  as  shown  in  Fig.  4.1-1;  data  for  4  March  1982 
1300  UTC  are  included. 

The  radiance  data  obtained  had  been  preprocessed  to  a 
'Level  1-B'  format.  Data  in  this  form  are  1)  quality  con¬ 
trolled,  2)  assembled  into  discrete  data  sets,  and  3)  un¬ 
calibrated  and  untransformed.  Calibration  coefficients  and 
earth  locations  are  included  in  the  data  set. 

The  data  contains  brightness  temperatures  which  must 
be  adjusted  using  band -correct ion  coefficients  and  converted 
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Figure  4.1-1  Coverage  of  Demonstration  Radiance  Data  Set 

to  radiances.  The  required  apparent  brightness  temperature, 
T* ,  is  given  by 

T  *  =  cTc  +  b  (4.1-1) 

where  Tc  is  the  band-corrected  brightness  temperature  and  b 
and  c  are  coefficients  dependent  on  the  sensor  and  channel. 
Using  T* ,  a  radiance  is  computed  using  the  Planck  Function 
(Eq.  2.1-2)  at  each  channel  frequency. 

With  a  grid  of  radiance  values  established,  radiance 
gradients  at  each  grid  point  are  computed  using  a  bicubic 
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spline  routine.  (Routines  from  the  IMSL  library  were  used  in 
this  effort.)  The  grid,  however,  is  aligned  with  the  satel¬ 
lite  scan  (s)  and  track  (t)  directions.  The  desired  gradients 
in  zonal  and  meridional  directions  are  given  respectively  by 


with 


dR 

dR 

dt 

dR 

ds 

dx 

dx 

ds 

dx 

dR 

dR 

dt 

dR 

ds 

dy 

dt 

dy 

ds 

dy 

dt 

dt 

dip 

dx 

dip 

dx 

dt 

dt 

dtp 

dy 

dtp 

dy 

dS 

ds 

dip 

dx 

dip 

dx 

ds 

ds 

dtp 

dy 

dy 

‘  V'. 

<P 

are 

earth 

longitude 

(4 .  1  -  2a) 

(4 . l-2b) 


(4 . l-3a) 


(4 . l-3b) 


(4 . 1 -3C) 


(4 . 1 -3d) 


spectively,  and,  under  the  assumption  of  a  spherical  earth, 

1 


dip 

dx 

dy 


R«  COS0 

J_ 

K 


(4 . l-4a) 

(4 . 1 -4b) 


where  R^  is  the  radius  of  the  earth  in  kilometers 


The  satel 

lite  track/scan  geometry  relative  to  earth  is  such  that  ana- 


lytic  expressions  for  the  derivatives  injiiij.',  >>i/c)(p,  <is/dip,  and  ds/f)<p 
do  not  exist.  However,  earth  latitude  and  longitude  for  every 
radiance  data  point  is  included  in  the  data  set.  Therefore 
these  derivatives  were  computed  numerically  using  finite  dif¬ 
ferences  . 

Finally,  to  compute  dR/ht  and  BR/ds,  the  radiance  data 
were  fit  to  a  surface  using  a  cubic  spline  routine.  The  de¬ 
rivatives  were  evaluated  from  this  surface  at  every  grid 
point . 

4.2  RAWINSONDE  DATA 

Rawinsonde  observation  (RAOB)  data  were  ordered  from 
the  National  Climatic  Data  Center  (NCDC)  to  provide  'truth' 
data  for  comparison  with  estimation  results.  Data  from  sev¬ 
eral  RAOB  sites  on  or  near  the  Iberian  Peninsula,  an  area  cov¬ 
ered  by  the  radiance  data,  were  ordered.  However,  for  4  March 
1982  1200  UTC  (the  synoptic  time  nearly  coincident  with  the 
collection  of  the  radiance  data)  upper  air  data  for  only  one 
station,  Gibraltar  (08495)  were  available.  This  finding  was 
both  disappointing  and  puzzling  since  Ref.  7  presents  maps  of 
thickness  values  at  various  pressure  intervals  for  these  Ibe¬ 
rian  s  ites . 

The  Gibraltar  u  and  v  wind  component  profiles  are 
shown  in  Fig.  4.2-1.  A  northwesterly  flow  from  1000  mb  to 
500  mb  is  indicated.  Above  this  layer,  the  flow  weakens,  and 
winds  become  light  and  variable  in  direction. 
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WMD  MAONTTUDE  (m/s) 

Figure  4.2-1  Gibraltar  Rawinsonde  Wind  Data 


4.3  TEMPERATURE  FIELD  DATA 

Paper  copies  of  850,  700,  500,  300,  200  and  100  mb 
National  Meteorological  Center  (NMC)  height  contour  and  tem¬ 
perature  analyses  were  ordered  from  NCDC  for  4  March  1982 
1200  UTC.  Wind  speed  analyses  replaced  temperature  analyses 
for  the  300  and  200  mb  levels.  Therefore,  these  temperature 
fields  were  subjectively  analyzed  over  the  area  of  interest. 
Temperatures  were  interpolated  to  a  one  degree  by  one  degree 
grid  for  display  and  gradient  calculation  purposes.  Fig¬ 
ures  4.3-1  through  4.3-4  show  the  850,  500,  200,  and  100  mb 

temperature  analyses. 
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The  850  mb  temperature  analysis  (Fig.  4.3-1)  indi¬ 
cates  a  strong  temperature  gradient  oriented  north-south  over 
the  Iberian  Peninsula.  The  isotherm  pattern  corresponds  to  a 
"dry"  cold  front  oriented  northeast  to  southwest.  The  500  mb 
isotherms  (Fig.  4.3-2)  indicate  a  strong  gradient  over  the 
northern  half  of  Spain;  the  gradient  weakens  southward.  At 
200  mb  (Fig.  4.3-3),  an  interesting  temperature  pattern  devel¬ 
ops.  A  cold  area  is  noted  near  Gibraltar  with  a  weak  tempera¬ 
ture  gradient  to  the  north  and  strong  gradient  over  Morocco 
and  Algeria.  Finally,  the  100  mb  isotherm  pattern 
(Fig.  4.3-4)  reveals  a  reversal  in  the  direction  of  the  gradi¬ 
ent  with  cold  air  to  the  south  and  warmer  air  to  the  north. 
This  pattern  is  expected  due  to  the  increasingly  higher  level 
of  the  tropopause  with  decreasing  latitude.  These  analyses 
along  w'ith  the  700  and  300  mb  analyses  (not  shown)  are  used  to 
compute  temperature  gradients  near  and  around  Gibraltar. 


5 .  BASELINE  SIMULATION  RESULTS 

For  purposes  of  algorithm  and  software  verification, 
simulation  results  have  been  generated  using  the  baseline  al¬ 
gorithm  (Section  3.1).  The  objective  here  is  to  analytical  1 v 
compute  radiance  gradients  from  a  simple,  prescribed  tempera¬ 
ture  field  and  the  satellite  weighting  function  models,  exer¬ 
cise  the  algorithm  with  these  simulated  data,  and  compare  the 
computed  thermal  wind  estimates  with  the  'true'  values  com¬ 
puted  directly  from  the  thermal  wind  equation.  The  simulation 
procedure  for  the  zonal  wind  component  is  detailed  below,  fol¬ 
lowed  by  simulation  results. 

For  simplicity,  but  with  no  loss  in  generality,  the 
temperature  field  is  prescribed  with  a  horizontal  gradient 
which  does  not  vary  with  altitude.  That  is, 

T(y,z)=(3y  +  7z  (5-1) 


where 


(3  =  =  CONSTANT 

ay 


(5-2) 


The  overbar  emphasizes 
zonal  thermal  wind  for 
Eq .  3.1-6  as 


the  constancy  of  the  derivative.  The 
the  layer  ,  z^ )  is  computed  from 


u 


T 


-R  ^ 

f  ay 


A  z 


(5-3) 


(i-J3 


where 


Az  =  z^-  z^  (5-4) 

This  represents  the  true  value  for  comparison  with  simula¬ 
tion  results . 

Radiance  gradients  are  computed  from  the  prescribed 
temperature  field  and  the  satellite  weighting  functions  (which 
depend  on  the  vertical  temperature  profile)  using  Eq.  3.1-4. 
The  computation  is  simplified  as  a  result  of  the  vertically 
constant  temperature  gradient.  At  each  relevant  wave  number, 

ILi  ^  I  K  .  dz  (5-5) 

dy  dy  J  ^ 

Since  the  are  represented  as  analytic  forms  (see  Sec¬ 
tion  2.1),  an  IMSL  routine  (DCADRE)  readily  computes  the  inte¬ 
grals  to  a  specified  accuracy.  The  computed  gradients  at 
selected  frequencies,  along  with  the  prescribed  temperature 
information,  are  the  inputs  to  the  baseline  algorithm. 

To  show  that  the  baseline  algorithm  will  produce  the 
correct  wind  estimate  in  the  absence  of  noise  and  numerical 
error  (due  to  using  the  trapezoid  rule),  let 


and  use  Eq.  5-5  to  get 


g  =  4^  L  (5-7) 

—  ay  - 

The  zonal  thermal  wind  estimate  is  given  by  (see  Eq .  3.1-15) 


-R  g"^  c 

(5-8) 

=  f 

where,  following  the  treatment  of  Section  3.1, 

SI  =  P-1  q 

(5-9) 

P  =  r  -t-  Of 

(5-10) 

q  =  aO 

(5-11) 

and,  for  this  case, 


aT  , 

—  L 
ay 

(5-12) 

^  A 

—  Az 

(5-13) 

ay 

Using  Eqs .  5-7,  5-8,  5-9  and  5-12,  the  thermal  wind  estimate 
for  vanishing  radiance  gradient  errors  becomes 


lim 

r-*o 


{uj} 


(5-14) 
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Applying  the  Sherman-Morr i son  formula  (Ref.  9)  to  Lq.  5-10, 
is  given  by 


-1  -  f - 

U  +  fr'oj 


(5-15) 


which  becomes 


lim  ,  lim  f-Rf,  O’^r'O 

r-.o<“T}  -  r-,()  (—('  ■  TTlfr^r-  "I 

=  lim  f^l' - a  ~|gTr-ig| 

r-0  \  f  1^1  +  e^r^ej^  -J 


-  R  -  R  ^»T  ^ 

-  «  = - Az  =  U  r 

f  f  ay 


(5-17) 


(5-18) 


(5-19) 


Thus  the  thermal  wind  estimate  ux  approaches  the  true  value  u^ 
as  the  radiance  gradient  error  covariance  matrix  approaches 
zero.  This  result  not  only  provides  theoretical  justification 
for  the  baseline  algorithm  but  also  enables  a  check  for  verti¬ 
cal  discretization  and  numerical  integration  errors. 
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Using  a  vertically  constant  meridional  temperature 
gradient  of  -0.0065  deg  K/km  and  estimating  the  zonal  thermal 
wind  for  the  layer  850-700  mb,  simulation  results  have  been 
generated  and  are  summarized  in  Table  5-1.  It  is  clearly 
shown  that  as  the  radiance  gradient  error  covariance  matrix, 
r,  approaches  zero,  the  estimates  approach  the  correct  value, 
as  dictated  by  the  theoretic  treatment  above.  Since  radiance 
gradient  errois  are  the  only  random  errors  presently  modeled 
in  the  baseline  algorithm,  the  estimate  root-mean-square-error 
also  vanishes  as  f  decreases. 

The  simulation  procedure  not  only  serves  to  verify 
the  algorithm  and  implementation,  but  also  provides  a  useful 
tool  for  sensitivity  studies.  For  example,  estimate  and  esti¬ 
mate  error  sensitivities  to  the  following  quantities  can  be 
determined 

•  Radiance  measurement  error  model  parameter 

•  Radiance  gradient  errors 

•  Satellite  weighting  functions 


TABLE  5-1 

BASELINE  ALGORITHM  SIMULATION  RESULTS 


ZONAL  WIND, 

(m/s) 

ROOT - MEAN -  SQUARE  -  ERROR 

(m/s) 

TRUE  VALUE 

3 . 5 

- 

BASELINE  RADIANCE 
GRADIENT  ERROR 
MODEL,  r 

2 . 5 

1 . 6 

r/10 

3 . 3 

0 . 7 

r/100 

3 . 4 

0 . 2 
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•  Temperature  errors 

•  Temperature  gradient  errors 

•  Channels  selected 

Further  work  in  this  area  will  assist  in  error  model  develop¬ 
ment,  particularly  for  the  satellite  weighting  functions  and 
temperature  and  temperature  gradient  errors. 

The  simulation  results  are  encouraging,  but  in  proc¬ 
essing  actual  radiance  data  the  wind  estimates  may  deviate 
from  observed  values  for  the  following  reasons 

•  nonzero  radiance  gradient  estimation  (i.e., 
calculation)  errors 

•  inexact  satellite  weighting  functions 

•  differences  between  true  winds  and 
geostrophic  winds 

•  observation  errors 

•  temperature  and  temperature  gradient  errors 

Demonstration  results  using  the  data  sets  described  in  Sec¬ 
tion  4,  among  others,  will  be  presented  in  the  final  report. 


6  . 


An  optimal  thermal  wind  estimation  procedure  directly 
utilizing  f requency - dependent  spatial  radiance  gradients  has 
been  described  herein.  The  procedure  is  optimal  with  respect 
to  the  models  developed  for  radiance  gradient  errors  and  sat¬ 
ellite  kernel  functions.  A  nominal  temperature  profile  and 
temperature  gradient  profile,  to  be  provided  operationally 
from  a  forecast,  are  key  elements  of  the  procedure.  This  in¬ 
terim  report  presents  both  a  baseline  formulation  and  an  ex¬ 
tension  of  that  approach  based  on  iteratively  updating  the 
temperature  data  along  with  the  thermal  wind  estimates. 

The  baseline  algorithm  and  software  was  verified  us¬ 
ing  simulated  radiance  gradients.  Simulation  results  were 
shown  to  be  consistent  with  baseline  algorithm  theoretical 
calculations.  Actual  HIRS/2  radiance  data,  described  herein, 
will  be  used  for  algorithm  demonstration  later  in  the  contract 
period,  and  results  presented  in  a  final  report. 
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APPENDIX  A-1 


The  objective  of  this  appendix  is  to  find  a  2  x 
Nx  *  Ny  matrix,  F.  ,  such  that  the  radiance  gradient  at  grid 
point  (j.k)  can  be  obtained  by  multiplying  the  vector  r  of  all 
radiance  measurements  over  the  grid  (as  defined  by  Eq.  2.2-2) 
by  F.  . .  This  form  is  useful  in  developing  the  radiance  gradi- 

J  >  ^ 

ent  error  mode  1 . 


Figure  A-l  illustrates  an  N  by  N  radiance  grid  with 

X  y 

horizontal  separation  distance  parameters  Ax  and  Ay.  The 
partial  derivative  of  radiance  r  at  point  (j.k)  with  respect 
to  X  is  calculated  using  r  values  along  the  j’^^  row  of  the 
grid;  the  partial  derivative  of  r  with  respect  to  y  is  calcu¬ 
lated  using  r  values  along  the  k*^^  column  of  the  grid.  The 
cubic  spline  routine  passes  a  set  of  cubics  through  the  radi¬ 
ance  values  of  either  row  or  column,  using  a  new  cubic  in  each 
subinterval.  In  order  for  both  the  slope  and  the  curvature  to 
be  the  same  for  the  pair  of  cubics  that  join  at  each  point, 
the  following  equations  must  be  satisfied  for  i  =  2,  3,  . . . , 

N^-1  (see  Ref.  8) 


-  Ax 


f 

V  c»x 


3rj,i  dr 

^  4  - + 


dx 


'n..u  "I 

dx  J 


and,  for  i  =  2,  3,  .  .  .  ,  -1: 
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Figure  A-l  Radiance  Grid 


Boundary  conditions  (known  as  "not-a-knot"  conditions)  are 
given  by; 


_Ax(4^  *2^]^  *  4r„. 


Av  I  2 — ~  +  4  — 

V  '^y 


.-,-2  +  4rj,.s-,_,  -  SFj.vj, 


5r,.k  -  4r2,k  -  r.i. 


^rv-i  .k  ,  3r 


^  3  ^ 
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Equations  A-1  through  A-6  permit  calculation  of  the  transfor¬ 
mation  matrix  F.  ,  .  Define  N  by  N  matrices  Q  and  P  as  the 

J  ,  K  XX  X  X 

coefficient  matrices  of  the  left  hand  sides  and  right  hand 
sides,  respectively  of  Eqs .  A-l,  A-3,  and  A-4.  Equivalently, 


For  instance,  i f  N.  =7,  then 

X 


2  4  0  0  0  0  0 
14  10  0  0  0 
0  14  10  0  0 
0  0  14  10  0 
0  0  0  14  1  0 
0  0  0  0  14  1 
0  0  0  0  0  4  2 


(A-8) 


(A-9) 


The  boundary  conditions  are  captured  in  the  first  and  last 
rows  of  and  ,  with  the  remaining  conditions  occurring  in- 
between.  Define  N  by  matrices  Q  and  P  similarly.  Then, 
if  ^  y  y 


The  final  step  consists  of  defining  an  (N^+N^)  by  (N^*  N^)  ma¬ 
trix  H  ,  such  that 
j  •  k 


It  can  be  easily  shown  that 


A- 1-4 


(A-14) 


■U 


where  denotes  the  unit  vector  in  (N  *  N  )  dimensional 

1  X  y 

space .  Hence 


J.k 


(A-15) 


is  the  desired  transformation  matrix  and 


Vfj.k  =  Fj.k  r 


(A-16) 


as  was  to  be  shown. 
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